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Abstract 

Path integral formulations for gauge theories must start from the canonical formulation in order to obtain 
the correct measure. A possible avenue to derive it is to start from the reduced phase space formulation. In this 
article we review this rather involved procedure in full generality. Moreover, we demonstrate that the reduced 
phase space path integral formulation formally agrees with the Dirac's operator constraint quantisation and, 
more specifically, with the Master constraint quantisation for first class constraints. For first class constraints 
with non trivial structure functions the equivalence can only be established by passing to Abelian(ised) con- 
straints which is always possible locally in phase space. Generically, the correct configuration space path integral 
measure deviates from the exponential of the Lagrangian action. The corrections are especially severe if the 
theory suffers from second class secondary constraints. In a companion paper we compute these corrections 
for the Hoist and Plebanski formulations of GR on which current spin foam models are based. 
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1 Introduction 



Path integrals for scalar Quantum Field Theories (QFT) on Minkowski space are supposed to compute the 
S - Matrix for the Hamiltonian in question. A standard Folklore says that heuristically one should simply 
consider all "paths" between some initial and final scalar field configuration on a spatial hypersurfaces 
labelled by moments of time ti,tf respectively and integrate over the exponential of (i times) the action 
with "Lebesgue measure" in order to obtain the evolution kernel. More specifically, let f2 be a (normalised) 
ground state (provided it exists) for the Hamiltonian H on a Hilbert space TL, let Q be the configuration 
space of spatial scalar field configurations with "configuration Lebesgue measure" dq, let Qt it t f = x te[u,t f ]Q 
be the set of paths, let [Dq] = Iltefo t/ ] dqt the "path Lebesgue measure", then 



e iS H,t f [q,q]/H 



< ty, j^-miu ^ >n= _ — — — (1.1) 

J Qhtf [Dq] n[q(t f )] Q[q(U)] e 1 *^/* 

and St it tf is the classical (Lorentzian) action integrated over the time interval [i»,t/]. Here by Lorentzian 
action we mean the Legendre transform 

S[q, q] := extr p { / ' dt[pq - H(p, q)]} (1.2) 

where p denotes the momentum conjugate to q and H the Hamiltonian. This "theorem" is wrong for several 
reasons. 

First of all, from the mathematical point of view, there is no Lebesgue measure on infinite dimensional 
spaces. Therefore one would like to consider [DQ] exp(iS/h) as a (complex) measure on Q but this does not 
work because the modulus of a complex measure [2] is supposed to be normalisable which is obviously not the 
case here. If the Hamiltonian is bounded from below, it is therefore much more promising to consider, instead 
of the unitary group M — * B(7i); t ^ exp(itH/h) the contraction semigroup M+ — > B(7i); t exp(—tH/h). 
Here B(Ti. denotes the algebra of bounded operators on H. Under these circumstances another Folklore 
theorem states that 



Sg.t, Ml/* 



< e-^f-^ A > H = _ SE (1.3) 

J Qt [Dq]n[q(t f )]n[q( ti )]e S W^ h 
where now Sf.^ denotes the "Euclidean" action, that is, the Legendre transfornfj] 

S E [q, q] := extr p { dt[ipq - H(p, q)]} (1.4) 

Ju 

While even under these circumstances the partition function 



Z:= [ [Dq]e- sB / h n[q(t f )]n[q( ti )] (1.5) 

Jq 

diverges, under fortunate circumstances it is possible to assign to e~ sE [Dq] jZ a well defined measure 
theoretic meaning on a proper a— algebra Q (with repect to which S E is usually not even measurable). 
Whenever (|1.3|) can be made rigorous, it is called the Feynman - Kac formula OH]. 

However, as it is well known [5j, also from the physical point of view, (jl.ip or (jl.3p are wrong in 
general. This is because the strict derivation (see e.g. [6]) of, say 111. 3b requires a skeletonisation of the 



1 Usually one obtains the Euclidian action by Wick rotation t — * it. However, we insist on this definition because it does not 
rely on an analytic structure of the fields in the time coordinate which is not justified anyway. Our definition is formally correct 
also in circumstances where the Hamiltonian is not only quadratic in the momenta with constant coefficients, see below. 
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time interval [ti,tf] and corresponding resolutions of the identity in terms of (generalised) position and 
momentum eigenvectors. That is to say, a priori one has to consider the complex hybrid action 

S C [q,p] := dt[ipq-H(p,q)] (1.6) 

which is integrated over both momentum and configuration coordinates. If H depends on p only quadratically 
with constant coefficients, then one can perform the Gaussian integral and up to an (infinite) constant 
which drops out in the fraction (|1.3j) one arrives at the Folklore result. However, in more general situations 
the result is different. For instance, the Hamiltonian could still be quadratic in p but with q dependent 
coefficients which leads to a nontrivial modification of the "measure" [Dq]. More generally, however, the 
Hamiltonian may not be quadratic or even analytic in p in which case an exact configuration space path 
integral representation is not available, only a saddle point approximation is available (plus the corresponding 
perturbative treatment of the non Gaussian corrections). Notice that the saddle point approximation and 
subsequent integrating out of the momentum variables reproduces (jl.4j) up to a non trivial measure factor. 

So far we have only considered scalar QFT on Minkowski space and even here we saw that the only 
correct derivation of the path integral proceeds via the Hamiltonian formulation, as stressed for instance 
in [5]. Additional technical and conceptual complications arise when we consider gauge theories and/or 
other background spacetime metrics. The simplest problem occurs for Yang - Mills type of gauge theories: 
Here the action is gauge invariant and if the measure is anomaly free (is gauge invariant as well) then one 
should divide by the (in general infinite) volume of the gauge group in order to give sense to both numerator 
and denominator in (II. 3D , If one considers QFT on non stationary background spacetimes then no natural 
Hamiltonian and vacuum exists [7] and the conceptual status of the path integral as a means to calculate 
scattering amplitudes becomes veiled. Even more veiled the situation becomes for totally constrained systems 
such as General Relativity on spatially compact four manifolds admitting globally hypebolic metrics when 
there is no true Hamiltonian at all. In this case certainly also the notion of a Wick rotation breaks down 
which on Minkowski space allows to reconstruct the Lorentzian Wightman functions from the Euclidian 
Schwinger functions via the Osterwalder - Schrader reconstruction theorem pQ. Parts of the reconstruction 
theorem, namely the construction of a Hamiltonian and a Hilbert space from a measure satisfying a natural 
background independent generalisation of the OS axioms, can be generalised to background independent 
theories [8j. 

It transpires that especially in the context of realistic physical theories, that is, General Relativity coupled 
to (standard) matter, it is neither clear what the heuristic Ansatz (jl.ip or (11. 3p computes nor whether it is 
the correct formula for what it is supposed to do. One possibility to deal with these problems is to try to 
solve the constraints classically and then to quantise the reduced phase space equipped with the (pull back 
of the) Dirac bracket [6j. This can be done in two ways. The first option is to impose suitable gauge fixing 
conditions in order to render the sytem totally second class and then to quantise the corresponding pull back 
of the Dirac bracket together with the induced reduced Hamiltonian. The second option is to determine 
explicitly a sufficient number of Dirac observables and to quantise the symplectic structure induced by the 
Dirac bracket. While for rare examples independent means exist to determine those gauge invariants, for 
most systems the only practical way to determine a sufficient number of Dirac observables is via a choice 
of gauge fixing. Namely, as we will review in the next section (see also e.g. the appendix in the second 
reference of |10|). there is a one to one correspondence between a choice of gauge fixing and a preferred set 
of gauge invariant functions which generate the full algebra of gauge invariant functions. In that sense the 
two methods, gauge fixing and this so called relational approach, are completely equivalent. The method 
is physically very interesting because it not only provides a suitable algebra of gauge invariant objects but 
also a gauge invariant Hamiltonian which drives the time evolution of those invariants. Here the question 
of equivalence between different choices of gauge fixing arises. As we will review in the next section, the 
preferred algebras of invariants that one obtains via different choices of gauge fixing are isomorphic. Of course 
they differ in their physical interpretation but as Poisson algebras they are isomorphic, the physical quantum 
kinematics is not affected by the choice of gauge fixing. The difference arises in the physical Hamiltonian, 
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that is, in the quantum dynamics. The explicit form of the physical Hamiultonian as a function of the 
invariant generators of the algebra of gauge invariant functions depends absolutely sensitively on the choice 
of gauge fixing and therefore even classically the evolution of the invariants will differ drastically from 
each other for different choices. For some choices the Hamiltonian may be explicitly time independent and 
leads to a conservative reduced system, for others it may not be. Even more crucial the choice of gauge 
fixing becomes in the quantum theory. Already for finite dimensional systems, depending on the choice of 
gauge fixing the physical Hamiltonian and other composite invariants built from the generators of the gauge 
invariant algebra may have discrete or continuous spectrum [9]! Notice that we here talk about composite 
invariants that have the interpretation of a given non invariant / measured in terms of another non invariant 
T (the so called clock) . If we change T to T' then its spectrum may switch from continuous to discrete or vice 
versa even though we talk about the same / and about the same Hilbert space representation!. In infinite 
dimensional situations the choice of gauge fixing has an even stronger influence for not only do we have 
to find a representation of the generators of the algebra of observables but in addition that algebra should 
support the physical Hamiltonian. One way to read Haag's theorem [II] is that Hamiltonians with different 
interaction terms cannot be implemented on the same Hilbert space. Thus generically different choices of 
gauge fixing will force us to choose different representations. For instance one may want to construct a cyclic 
representation built from application of the generators to a vacuum (ground state of the Hamiltonian). That 
vacuum of course depends on the Hamiltonian and even for free field theories those cyclic representations 
are typically unitarily inequivalent. In case that the physical Hamiltonian is explicitly time dependent, one 
is in addition confronted with the usual problem of QFT on curved spacetimes, namely that one has to 
decide at which point of time one wants to select a vacuum vector. 

All of this certainly strongly affects the resulting reduced phase space path integral because it is based on 
the selected Hilbert space representation and the transition amplitudes between physical states do depend 
on the physical Hamiltonian. For any such choice, the path integral then does have the interpretation of 
(jl.ip or (|1.3j) in terms of the reduced Hamiltonian. 

An additional complication that we have not mentioned yet is the case of a system with second class 
constraints. Such a system is to be canonically quantised with respect to the Dirac bracket rather than the 
Poisson bracket. Typically the Dirac bracket destroys the canonically conjugacy of the global coordinates of 
the phase space that one started from. Since to find representations of such complicated Poisson algebras is 
usually prohibitively difficult one is forced to switch to local Darboux coordinates (by means of a canonical 
transformation with respect to the original Poisson bracket) which is always possible locally [6|. Such 
coordinates may be very difficult to find in practice. Assuming this to have done nevertheless, one can then 
construct the reduced phase space using a choice of gauge fixing as already described above and after having 
chosen a Hilbert space representation subordinate to that gauge fixing, the transition amplitudes in terms 
of the induced physical Hamiltonian. 

From here on then mostly one proceeds rather formally. One assumes that one can choose a Schrodinger 
representation based on the reeduced Darboux configuration space. By using well known skeletonisation 
techniques one then basically writes the transition amplitude between initial and final states \E'j,\I'/, as a 
path integral over the reduced Darboux phase space, replacing the reduced Hamiltonian operator by its 
classical function which results in the exponential of the reduced Hamiltonian Darboux action. In order 
to make contact with (jl.ip one wants to rewrite this path integral as a path integral over the unreduced, 
original configuration space and in terms of the original Lagrangian. As is well known, this can be formally 
done and we will review this rather involved procedure in section [3l Basically one first extends the reduced 
Darboux phase space to the unreduced Darboux phase space thereby introducing 5 distributions of the 
constraints and the gauge fixing condition as well as measure factors which cancel the Jacobian that arises 
when solving the 5 distributions. One then observes that, in presence of the 5 distributions the reduced 
Hamiltonian action can be written as the unreduced symplectic potential, in terms of the unreduced Darboux 
coordinates. Interestingly, the measure factors and the 5 distributions combine in just the right way as to 
make the resulting expression independent of the gauge fixing condition when considered as a measure on 
gauge invariant functions. This is similar to the Fadeev - Popov theorem [6] and we will review this result 
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in section This seems to be in contradiction to what we have said above about the dependence of the 
transition amplitudes on the gauge fixing condition. The resolution is that at this point the integral is not 
over gauge invariant functions, it is an integral over ^fi^ff which are functions at initial and final points of 
time of the reduced Darboux coordinates which are not gauge invariant. More generally, in applications to 
scattering theory, we may be interested also in n-point functions so that the path integral is over functions of 
the reduced Darboux coordinates also at intermediate times (in fact we will use the method of a generating 
functional so that there is a dependence on the reduced Daroux coordinates at all times). One may, 
in the presence of the 5 distributions, extend the non gauge invariant, reduced Darboux coordinates to 
gauge invariant functions which use the chosen gauge fixing condition. However, these extended functions 
now display a complicated dependence on all unreduced Darboux coordiantes which makes this extension 
practically useless. Even if one did perform the extension, while one can now change the dependence on the 
gauge fixing condition in the measure, one cannot get rid of it in the gauge invariantly extended functional. 
In any case, one next performs the canonical transformation that leads from the Darboux coordinates back to 
the original canonical coordinates which does not affect the symplectic potential and the Liouvile measure 
but it affects the initial and final states. Then one exponentiates the 5 functions and, by the technique 
introduced in [12] . gets rid of the secondary second class constraints which leads to further changes in the 
measure. Finally, one integrates out the momenta. This is only possible if the reduced Darboux configuration 
coordinates, as functions of the original canonical coordinates, do not depend on the original momenta and 
if they do not leads in general to further changes in the measure while now the exponential of the covariant 
Lagrangian action appears. 

The point of mentioning these in principle well known facts is twofold. The first is that we wish to stress 
that even if all the assumptions that we have listed can be verified, the correct Langragian configuration 
space measure may differ drastically from the naive one in These deviations depend crucially on the 

dynamical content of the theory and cannot be discarded. The second point that we want to make is the 
dependence of the transition amplitudes on the chosen gauge fixing. This dependence is at first astonishing 
because one is used from Yang - Mills theory that the path integral does not depend on the gauge fixing 
and it even sounds dangerous because it seems as if this dependence implies that gauge invariance is broken. 
However, this is not the case: The dependence on the gauge fixing is physically correct. The reason is that 
in generally covariant systems the dynamics mixes with gauge invariance. In Yang - Mills theory this is 
not the case, there one has a gauge invariant Hamiltonian at one's disposal which is not generated by a 
gauge fixing condition, it is simply there without further input. Gauge invariant functions in Yang - Mills 
theory can also be easily constructed without ever mentioning any gauge fixing, for instance Wilson loops or 
flux tubes between quarks. The gauge fixing condition comes in only when cancelling an otherwise infinite 
constant. This introduces a gauge fixing 5 distribution and a Fadeev - Popov determinant into the measure 
whose combination is independent of the gauge fixing by construction, similar as in our discussion above. In 
contrast, in generally covariant systems a gauge fixing condition can be seen as defining a preferred algebra 
of observables and a preferred dynamics thereof. Gauge invariance is not at all broken, the dynamical system 
consisting of reduced Darboux phase space and reduced Hamiltonian as defined by a gauge fixing is in one 
to one correspondence with a dynamical system consting of a preferred algebra of Dirac observables and 
a gauge invariant physical Hamiltonian defined via the same gauge fixing (now interpreted as a choice of 
clocks). The two descriptions are equivalent. The gauge fixing dependence comes in because one needs a 
gauge fixing in order to arrive at the very notion of a dynamics, or in other words, at the very notion of 
an observer. This observer dependence of the classical and quantum theory has already been stressed in 
[13j and will be discussed in more detail in |14j . Let us stress again, as we have already said, that similar 
as in Yang - Mills theories the gauge fixing dependence of the measure disappears when we restrict it (as 

2 A special situation arises if one considers gauge transformations that tend to the identity in the infinite past and future 
and that the only non gauge invariant functions in the path integral are located at the infinite past and future. This is not 
the case for the n-point functions or the generating functional but for the rigging kernel between two kinemtical states. Now 
the dependence on the gauge fixing formally disappears from the path integral, of course modulo the representation theoretic 
caveats that we have mentioned. 
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a linear functional) to gauge invariant functions. However, the choice of those gauge invariant functions 
themselves and the corresponding physical Hamiltonian, in other words the physical interpretation of the 
theory, induced by a choice of gauge fixing (clock) is what makes the description gauge choice dependent. In 
contrast, in Yang - Mills theories such a choice of clocks is not necessary in order to arrive at useful gauge 
invariant functions. In principle, the generators of the algebra of gauge invariant functions for one choice of 
gauge fixing can be written as complicated functions of the generators for any other choice. However, this 
involves an infinite series of commutator functions about whose convergence nothing is known and which 
therefore is practically useless if not mathematically ill -defined. 

In this paper we want to illustrate the complications sketched above for a general theory which will be 
the first result of this paper. While certainly bits and pieces of our description appear in various places 
in the literature, we hope that assembling them in the form presented here may add a certain amount of 
clarity to the question how reduced phase space and path integral quantisation fit together. 

The second result of this paper will be to sketch how the path integral is related to Dirac's operator 
constraint quantisation [15] and a particular incarnation of it, the so called Master Constraint Programme 
|16j for first class systems. As already mentioned, the reduced phase space rarely admits a global Darboux 
coordinate system and hence a quantisation of the unreduced phase space is much simpler. The price to pay 
is that one has to impose them as non - anomalous operators on that Hilbert space in order to compute the 
physical Hilbert space. There are certain heuristic group averaging methods [T7] available in the literature 
which, as the name suggests, apply when the constraints form a Lie algebra. If they do not (structure 
functions), then not only are the constraints difficult to define without anomalies because of factor ordering 
difficulties but also group averaging is not applicable. It is for that reason that the Master Constraint Pro- 
gramme (MCP) was introduced. In the MCP, all constraints are enconded into one single Master constraint. 
The Master Constraint is a classically equivalent platform and is automatically free of anomalies so that 
group averaging (or direct integral decomposition) methods apply. 

The central ingredient of the group averaging method is a (generalised) "projector" (or rigging map) from 
the kinematical Hilbert space into the physical one, equipped with an associated physical inner product. 
It can be expressed in terms of a path integral which in case of a true Lie algebra is readily recognised as 
(jl.ip or (jl.3p respectively. In case of the Master constraint that can also be established, however, the proof 
is somewhat more involved. Not surprisingly, the key to the understanding of how all of these methods fit 
together is how the reduced phase space description arises from the constraints and a suitable gauge fixing 
condition which in turn allows for a local Abelianisation of the constraints. It may seem astonishing that the 
gauge fixing condition enters the interpretation of the physical Hilbert space in such a prominent way. The 
reason for why that happens is that the physical Hilbert space can be considered as the closure of the set 
of vectors that one obtains by applying the algebra of gauge invariant observables to a cyclic physical state. 
However, the construction of that algebra and the interpretation of its elements is faciliated by considering 
the gauge invariant extension of the kinematical algebra as induced by a gauge fixing condition. In other 
words, while in the operator constraint method one only deals with manifestly gauge invariant objects, 
their interpretation again relies on a gauge fixing condition or equivalently on a choice of rods and clocks. 
Different such choices result in the same algebra but its generators (elementary observables) differ for each 
choice. 

To summarise: 

The correct path integral formula and its interpretation can only be obtained by following the Hamiltonian 
path, otherwise one misses important corrections to the measure. In the context of spin foam models [18] for 
Loop Quantum Gravity [19] this has been pointed already in [20] (see also [21], [22]). The corrections to the 
measure are not manifestly covariant as first indicated in [23] but seem to be required in order to maintain 
at least some form of spacetime covariance as claimed in [24 j. They should therefore be taken seriously in 
any realistic spin foam model for General Relativity. Work is now in progress which tries to implement these 
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corrections. See also [22] where the covariance of the path integral with respect to the Bergmann - Komar 
group is discussed. 



This article is organised as follows. 



In section two we review Dirac's analysis of gauge sytems and the relation between gauge fixing and a 
gauge invariant description. 

In section three we derive the path integral from the quantisation of the reduced phase space based on 
the natural generators of the algebra of invariants defined by a choice of gauge fixing conditions. 

In section four we derive the same path integral starting from the unreduced (with respect to the first class 
constraints) phase space and implement the rigging map. In case that the first class algebra of constraints 
involves structure functions, using the rigging map technique requires to pass to new constraints that form 
an algebra. This is always (locally) possible because one can always (locally) Abelianise constraints. 

In section five we use as an alternative route the MCP and show that again one arrives at the same 
path integral. This is to be expected because both constraint rigging and Master constraint rigging should 
provide the generalised projector on physical states, however, the technical mechanism by which this works 
is somewhat involved. 

Finally in section six we summarise and conclude. 



2 Classical Preliminaries: Gauge Fixing Versus Gauge Invariant Formu- 
lation 

In an attempt to make this article self - contained we start with the classical theory. We will neeed the 
corresponding notation anyway for the path integral formulation. First we summarise the main ingredients 
of Dirac's algorithm. Then we display the relation between the reduced phase space of gauge invariant 
observables and the pull back phase space as induced by a gauge fixing. As we will see, the two formulations 
are equivalent for suitable choices of gauge fixing. 

2.1 Brief Review of Dirac's Algorithm 

We consider a theory with Lagrangian L(q a (t),v a (t)) and corresponding action 

S = [ dtL(q a (t),v a (t)) (2.1) 
JR 

Here the index a takes values in a general set which may comprise discrete and/or continuous labels. We 
are interested in a theory with gauge symmetries so that the Lagrangian will be singular, that is, we cannot 
solve all the velocities v a = q a for the canonical momenta 

* = £ < 2 - 2 > 

By solving a maximal number of velocities v a (whose number is equal to the rank of the matrix d 2 /dv a dv b ) in 
terms of the momenta p a and the remaining velocities v 1 , that is, v a = u a (q a ,p a ;v l ) such that (v a ) = (v a ,v l ) 
(i.e. the indices a and i take values in index sets that partition the index set associated with a) we obtain 
the primary constraints 

dL 

Ci = pi — [-Q^-} v <x=u a (2-3) 
which does not depend on the v 1 by assumption of the maximality of the the v a . The canonical Hamiltonian 

H c = [v a p a - L{q, v)] va = u « (2.4) 
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always has the structure [EJ [H5] 

H c = H^q,p)+v i C i (q,p) 



(2.5) 



that is, it is an affine function of the v l . 

The further analysis of the system is now governed by Dirac's algorithm |15j : 
One requires that the constraints are preserved by the Hamiltonian flow of H c . Whenever {H c ,Ci} = is 
not satisfied on the constraint surface, there are two possibilities: Either 1. {H c ,Ci} does not involve the 
velocities v l or 2. it does. In the first case we must add {H c , C{\ to the list of constraints in the second we 
solve all the equations of type 2 for some of the velocities v l (assuming that the system of equations is not 
overdetermined) . Iterating like this, one ends up, in general, with further constraints Cj, which are called 
secondary constraints, and the velocities are restricted to be of the form v % = VQ(q,p) + X m v t m (q,p). Here 
v l = Vq solves {H c , Cj} = {H c , Cj} = for all j, J on the constraint surface {Cj = Cj = OVj, J} and (v l m ) m 
is a maximal linearly independent set of solutions of the system v l {Ci, Cj} = v l {Ci, Cj} = for all j, J (on 
the constraint surface). The coefficients X m are free and phase space independent. 

It follows that the 

F m := via (2.6) 

are first class constraints, i.e. they weakly (i.e. on the constriant surface) Poisson commute with all 
constraints. By taking linear combinations of the constraints Ct, Ci (with phase space dependent coefficients) 
we isolate a maximal number of first class constraints. The constraints F m are called primary first class 
constraints, the additional ones Fm are called secondary first class constraints. The remaining constraints 
among the set (Cj,C/) which are linearly independent of the set (F^) := (F m ,FM) are called second class 
constraints and are denoted by (<Se)- 

The canonical Hamiltonian can now be written 

H c = H% + X m F m (2.7) 

where 

HS = H' + vlC, (2.8) 

is of first class by construction. It may therefore be an affine function of the F^ (with phase space dependent 
coefficients) 

H{ = H + f^F^ (2.9) 

The piece Hq is referred to as the true Hamiltonian because it is not constrained to vanish. In totally 
constrained systems such as General Relativity it vanishes identically, that is, the canonical Hamiltonian is 
a linear combination of first class constraints. Notice that only the primary first class constraints appear in 
the canonical Hamiltonian with arbitrary coefficients X m and so one would associate gauge invariance only 
with respect to them. However, this is in general inconsistent because the Poisson algebra of primary first 
class constraints generically does not close, only the full set of first class constraints always does. In other 
words, since the Poisson bracket between two first class functions is first class and if O is weakly invariant 
under the F m then also 2{F[ m , {F n i , O}} = {{F m , F n },0} should be weakly zero. Furthermore, the time 
evolution of O with respect to H c should be gauge invariant which gives 

{F m , {H c , O}} = -{H c , {O, F m }} - {O, {F m , H c }} « {{F m ,Hg}, O} « (2.10) 

Hence gauge invariant quantities should be those that weakly Poisson commute with the minimal subset of 
first class constraints generated by the Poisson brackets between the {H c , F m } (and higher order brackets 
with H c ) and between the first class primary constraints among each other. For most systems of physical 
interest this exhausts all first class constraints and we will assume this to be the case here. In that situation 
the piece H" of the Hamiltonian in (12. 9p will therefore generically contain all secondary first class constraints 
as well, that is, the corresponding phase space functions will be non - vanishing. See [U [27] [28] for a 
discussion when this so called Dirac conjecture can be proved. 
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As far as the second class constraints are concerned, they are not associated with any gauge freedom. 
It is in fact inconsistent in general to require an observable to satisfy {Ss, O} as an application of the 
Jacobi identity reveals. This means that observables are not first class functions, they only have to weakly 
Poisson commute with the first class constraints, not with the second class constraints. Hence, to solve 
the second class constraints we simply have to restrict ourselves to the corresponding constraint surface. 
In other words, once we have computed the functions on phase space which have weakly vanishing Poisson 
brackets with all first class constraints, we should restrict them to the constraint surface defined by the 
second class constraints only. The induced symplectic structure between such observables /, /' on the total 
constraint surface is simply the pull back of the symplectic structure on the unconstrained phase space by the 
embedding of the constraint surface defined by the second class constraints into the full phase space. More 
precisely, let V denote the unconstrained phase space and V := {m 6 V\ S^(m) = VS} the constraint 
surface defined by the second class constraints. Consider the corresponding embedding J : V — * V- If 
denotes the symplectic structure on V then 0~ := J*Q, denotes the pull - back symplectic structure on V . 
This is again a symplectic structure because it is closed and non degenerate which follows from the fact that 
the matrix 

A ss , :={5 S ,5 E /} (2.11) 
is non degenerate. The corresponding Poisson bracket is given by 

{J*f,rf}- = r{f,f}* (2.i2) 

where 

{/,/'}* :={/,/'} -{/,5 S } (A- 1 )^' {Sv,f'} (2.13) 

denotes the Dirac bracket on the full phase space. We will prove this for the convenience of the reader in 
the next section. 

The reduced phase space is defined by the Poisson algebra of gauge invariant observables, which are 
not weakly vanishing (i.e. which are not linear combinations of first class constraints) equipped with the 
pull-back of the Dirac bracket to the constraint surface defined by the second class constraints. Notice 
that the Dirac bracket generically changes the symplectic structure for the observables as compared to the 
Poisson bracket. However, it does not change the equations of motion defined by the canonical Hamiltonian 
as the Dirac bracket and the Poisson bracket between two functions coincide whenever at least one of them 
is a first class function. 

2.2 Reduced Phase Space and Gauge Fixing 

In principle the description of the previous subsection entails the complete information about the dynamics 
and the physical (gauge invariant) content of the theory. However, it does not provide an explicit description 
of the observables. Moreover, in totally constrained systems the equations of motion for the observables 
with respect to the canonical Hamiltonian are trivial which means that some important ingredient is missing 
in that case: A non vanishing physical Hamiltonian which drives the time evolution of the observables. In 
this section we give an explicit construction of the reduced phase space, provide a physical Hamiltonian and 
display the relation of our framework to gauge fixing. 

We saw that we eventually obtain a constrained Hamiltonian system with first class constraints and 
second class constraints S% on a phase space with canonical pairs (q a ,p a ), a = 1, ..,n; m < n with respect to 
the original Poisson bracket. As shown in [27], there always exists a local canonical transformation (with re- 
spect to the Poisson bracket) from the canonical pairs (q a ,p a ) to canonical pairs (Q A , Pa), (4^ '>7T/t)> ( x<7 ' iVo) 
such that 

5 S = ^ z s = (2.14) 

where the index a takes half the range of that of S and where we denoted either x a or y a by for some 
value of S. It is then clear that at least weakly the Dirac bracket and the Poisson bracket coincide on the 
(Q, P, 7r, 4>) and that z has zero Dirac bracket with anything. 
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Next, if there is a true, gauge invariant Hamiltonian Hq in (|2.9p (not constrained to vanish) enlarge 
the phase space by an additional canonical pair (q ,Po) and additional first class constraint Fq = po + H. 
The reduced phase space and dynamics of the enlarged system is equivalent to the original one, hence we 
consider without loss of generality a system with no true Hamiltonian (totally constrained system). The 
canonical Hamiltonian of the system is then a linear combination of the first class constraints 

H c = p^F^ (2.15) 

Here we have set p M = f for secondary first class constraints and p m = f m + A m for primary first class 
constraints where is defined in (|2.9[) . 

A gauge fixing is defined by a set of gauge fixing functions with the property that the matrix with 
entries M^ u := {C^,G V } has everywhere (on the unconstrained phase space) non vanishing determinant^. 
Notice that we allow for gauge fixing conditions that display an explicit time dependence. The conservation 
in time of the gauge fixing conditions 



TO r\ 

jG» = q-G^ + {H c , G^} « -G„ + p v M Vil (2.16) 



uniquely fixes the "Lagrange multipliers" to be the following phase space dependent functions 

p " = {M ~' r & (2 - 17) 

At this point one may be puzzeled by the following issue: The functions p^ already depend on the phase 
space through For the p m we can always solve (|2.17p for the free function A m . But for the p M the 
solution (|2.17p leads to a consistency condition on the already imposed gauge fixing conditions, in other 
words we should impose independent gauge fixing conditions only for the primary first class constraints^. 
This is indeed true as far as fixing the free coefficients in the canonical Hamiltonian is concerned. However, 
in view of the fact that all first class constraints generate gauge transformations, one has to eventually 
reduce with respect to all their gauge motions. Therefore it is mathematically and physically equivalent 
and mathematically much more convenient to regard all p^ as free parameters, that is, to drop the phase 
space dependence of the f M . Hence to fix the gauge we need gauge fixing conditions for all first class 
constraints. We will see explicitly in the path integral formulation that one is forced to this point of view 
and that nevertheless one can restore the phase space dependence of the f M when eventually reducing the 
path integral as one over configuration space (rtaher than the phase space). 

By construction of the Dirac bracket, we can simply ignore the variables z for what follows and set them 
equal to zero where ever they occur. In terms of the remaining canonical pairs we can solve Fn = G^ = 
for 

F IM = TT ll + h tl (Q,P) = 0, G„ = ^-t»(Q,P) = (2.18) 



3 Ideally, the gauge G M = should define a unique point in each gauge orbit. 

4 As an example, in General Relativity the primary constraints demand that the momenta conjugate to lapse and shift vanish, 
the secondary constraints are the spatial diffeomorphism and Hamiltonian constraints respectively. All constraints are first class 
and the canonical Hamiltonian is a linear combination of all of them, in particular lapse and shift play the role of the f M for the 
Hamiltonian and spatial diffeomorphism constraint respectively. A consistent gauge fixing would now be to first prescribe four 
functions purely built from the intrinsic metric and their conjugate momenta (independent of lapse and shift). Such conditions 
have vanishing Poisson brackets with respect to the primary constraints. Therefore equation (|2.17|l can be computed and 
presribes lapse and shift as a function of intrinsic metric and conjugate momentum alone. The remaining four gauge fixing 
conditions for the velocities (Lagrange multipliers) of lapse and shift which are the coefficients of the primary constraints are 
now that they are the time derivatives (Poisson brackets with the canonical Hamiltonian) of the already prescribed functions 
for lapse and shift. These conditions are then consistent with the equations of motion, i.e. that the Lagrange multipliers are 
the time derivatives of lapse and shift. The corresponding matrix {J)j,G„} in this case is block diagonal. We could also have 
prescribed lapse and shift in the first place as functions of intrinsic metric and conjugate momentum and then would have to 
find four additional gauge fixing conditions on those variables whose equations of motion lead to the prescribed values of lapse 
and shift. 
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for certain functions h, r which generically will be explicitly time dependent. The variables 4>, tt are called 
the gauge degrees of freedom and Q, P are called the true degrees of freedom (although typically neither of 
them is gauge invariant). 

The reduced Hamiltonian H red (Q,P), if it exists, is supposed to generate the same equations of motion 
for Q, P as the canonical Hamiltonian does, when the constraints and the gauge fixing conditions are satisfied 
and the Lagrange multipliers assume their fixed values (|2.17p . that is, 

{#red, /} = {#can, f}F=G=p-po=0 = [Po{Ffi, f}]F=G=p-po=0 (2.19) 

for any function / = f(Q,P). For general gauge fixing functions the reduced Hamiltonian will not exist, 
the system of PDE's to which (|2,19p is equivalent to, will not be integrable. 

However, a so called coordinate gauge fixing condition = (p^ — with independent of the phase 
space always leads to a reduced Hamiltonian as follows: We can always (locally) write the constraints in the 
form (at least weakly) 

Fp = AV(7r„ + h!M Q, P)) =■ Mpv K (2-20) 

where h^(Q,P) = h' ((f) = r,Q,P). Notice that the locally equivalent constraints F'^ are actually Abelian 
by a general argument [6]. Then, noticing that M^ v « {F^,G U }, (|2.19p becomes 

{H rcd , /} = [p^M^lK, f}] F =G= P -p =o = [pp{hp, f}] G =o = {Pih'^, /} (2.21) 

with hn = h^(<p = t, Q, P) and we used that / only depends on Q, P. This displays the reduced Hamiltonian 

as 

H Icd (Q, P; t) = r"(i)M0 = Q> P )) ( 2 - 22 ) 

It will be explicitly time dependent unless rj is time independent and is independent of cf>, that is, unless 
those constraints can be deparametrised for which 7^ 0. Hence, deparametrisation is crucial for having a 
conserved, reduced Hamiltonian system. 



On the other hand, let us consider the gauge invariant point of view. Following the general framework 
[291 l30l EH [32l EH E2] it is possible to construct a gauge invariant extension of any gauge variant function 
f(Q,P) off the gauge section <j> = r by the following formula 

O f (r) = [expWXj ■ ftp^r-t (2.23) 

where we have denoted the Abelian Hamiltonian vector fields by X^ := {ir^ + hn, •}• It is easy to check 
that {Of(r),FfjL} ~ 0. Consider a one parameter family of flows t *— ► t^(£) then with Of it) := Of(r(t)) we 
find 

(j rhi m*. 

-O f (t) = f»(t) £ X.X^.X^ ■ f (2.24) 

n=0 

On the other hand, consider H^(t) := O^ (r(t)), then [32] 

{#„(*), 0,(t)} = =0 {h ^f Y {T{t)) = {K j } (r(t)) = O x ^.f(T{t)) 

°^ mi RUn 

= X.X^.X^.f (2.25) 

n=0 

Here the bracket {.,.}* denotes the Dirac bracket associated with the second class system (F^^G^). In the 
second step we used that neither nor / depend on ir u , in the third we used that / does not depend on <f>" 
and in the last we used the commutativity of the X u . Thus the physical Hamiltonian that drives the time 
evolution of the observables is simply given by 

H(t) := t^t)h'(r(t), Q (t), P (t)) (2.26) 
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where we used that (|2.23|) is a Poisson automorphism [32], that is, 



{O f (r),O g (T)} = {M (r) (2.27) 

Here {f,g}* = {f,g} for functions of Q,P only was exploited. This is exactly the same as (j2.21|) under the 
identification / *-* Oj(0). Hence we have shown that for suitable gauge fixings the reduced and the gauge 
invariant frameworks are equivalent. Notice that it was crucial in the derivation that (^,7r^) and (Q , Pa) 
are two sets of canonical pairs. If that would not be the case, then it would be unclear whether the time 
evolution of the observables has a canonical generator. 

The power of a manifestly gauge invariant framework lies therefore not in the gauge invariance itself. Rather, 
it relies on whether the gauge fixing can be achieved globally, whether it can be phrased in terms of separate 
canonical pairs, whether the observer clocks cj^ are such that reduced Hamiltonian system is conserved and 
whether they do display the time evolution of observables as viewed by a realistic observer. See [101 PH] for 
a discussion of this point. 

Our description sketched above shows that a useful, manifestly gauge invariant formulation implicitly 
also relies on a system of gauge fixing conditions. Namely, the gauge fixing conditions equip us first with 
a physical interpretation of the Dirac observables and second with a physical time evolution: 0/(r) has 
the meaning of that relational observable, i.e. gauge invariant quantity, which in the gauge (j) = r takes 
the valu^l /. Its time evolution in terms of r is driven by a gauge invariant Hamiltonian which arises by 
solving the constraints F = in terms of the momenta ir conjugate to <f>. The gauge fixing condition thus 
prominently finds its way into the very interpretation of the physical (reduced) phase space. If we would 
choose different clocks <f>' then different observables 0'^{t) would result. Due to (|2.27p . the algebra of the 
0/(t) among each other and of the 0'^{t) among each other respectively are isomorphic provided that / 
and /' only depend on the respective true degrees of freedom. In particular, both the (Oqa(0), Op A (0)) 
and (O'q IA (0),O' p , (0)) respectively provide a (local) system of coordinates on the reduced phase space and 

therefore one can translate between the twcH. However, their physical interpretation and physical time 
evolution is entirely different. This crucial fact will also be reflected in the interpretation of the path 
integral. 

3 Reduced Phase Space Path Integral 

This section is subdivided into three parts. In the first we make some general remarks about scattering 
theory in ordinary QFT which is closely related to the path integral and how this applies to our case. In 
the second we formally derive the reduced phase space integral as the generating functional of n— point 
functions. This path integral is an integral over the reduced phase space. In the third section we unfold this 
path integral and integrate over the unconstrained phase space whereby proper gauge fixing conditions and 
constraints have to be imposed. 

5 One maybe tempted to run the following contradictory arguement: 0/(r) obviously coincides with / in the gauge <f> = r. 
Since it is also gauge invariant and since any other gauge can be reached from (j> — r one may think that it takes the value / in 
any other gauge, say (f> = t' which is obviously not the case by inspection. The catch is that in order to reach the gauge <f> = t 
from cj> = t one must apply a gauge transformation to 0/(r) which maps <j> to (f> + r' — r — (f> + St and / to its corresponding 
image f + Sf under this gauge transformation. By gauge invariance we obtain 0/(r) = 0/+a/(r + St). Hence in the gauge 
<j> = t' the observable takes the value f + Sf and not /. This is not in contradiction with gauge invariance because / + Sf and 
/ are evaluated at different points on the the same gauge orbit just in the right way as to give the same numerical value. 

B Notice that when choosing different clock variables <t>',n' we also have to choose different true degrees of freedom Q',P'. 
The algebra of the Oq(0), O' p (0) is not isomorphic to the one of Oq(0), Op(0), rather we have {O' P (0), O' Q (0)} = 0' {PQ} „ (0) 
where {.,.}'* denotes the Dirac bracket associated with F, <j>' . Thus, while Op(0), Oq(0) and O'p, (0), O'q, (0) are conjugate pairs, 
O' P (0), O'q (0) are not. 
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3.1 Remarks about Scattering Theory 

The central object of interest in QFT is the scattering matrix. Rigorous scattering theory is in fact a difficult 
subject even in ordinary QFT on Minkowski space. First of all, there is a notion of a free and interacting field 
<p and (j) respectively which evolve according to the free and interacting Hamiltonian Hq and H respectively. 
Here free means that Hq does not contain any self - interaction. The physical assumption is that in the 
far future tf — > oo and far past t{ — ► — oo any outgoing and ingoing particles respectively do not interact. 
This is, of course, not really true. However, using the methods of local quantum physics, assuming that the 
theory has a mass gap0 one can prove that the vacuum correlators of the asymptotic fields reduce to those 
of the free field, where vacuum really means the interacting vacuum. 

This means that the asymptotic fields generate from the interacting vacuum a Fock space 7i± which in 
general could be a proper subspace of TL. These states can be thought of as the rigorous substitutes for 
the states generated by the non existing asymptotic free field from the free vacuum. This is the famous 
framework of Haag and Ruelle, see [TT] and references therein. The rigorous S - matrix is then defined 
by the scalar product between these asymptotic Fock states which one interprets as vector states in the 
Heisenberg picture under the free dynamics. The rigorous relation between the S - matrix elements and 
the time ordered n— point functions is then provided by the famous LSZ formula [11] which rests on the 
assumption of asymptotic completeness, that is, TL = TL±. 

In most textboks on QFT, the Haag - Ruelle theory is barely mentioned. Rather, one somehow postulates 
that the free and interacting fields can be implemented on the same Hilbert space. Then one can consider 
Fock states ip± = i/j± created from the free vacuum by free field creation operators which one considers 
as vectors in the Heisenberg picture in the infinite past and future respectively because they are time 
independent. To obtain the corresponding Schrodinger picture states in that limit one just has to evolve via 
the free Hamiltonian to obtain ip±(t) = e lH ot/^^H ag ^ _^ -)- 00 _ ^0 map a Schrodinger state from the far past 
to the far future one should however use the interacting Hamiltonian H and get the evolved Schrodinger 
picture vector state e lH ^ t f~ ti ^ ti i[)^.(ti). Thus the scattered Heisenberg picture state would be given by 

Mattered =lim S(t + ,t-)^ (3.1) 

where 

S(t f ,ti) := V{t f ) ] V(U) V(t) = exp{-itH/h)exp{itH /h) (3.2) 

The formally unitary operators V(t) in principle map the evolving free Heisenberg field operators to the 
evolving Heisenberg field operators. Using the differential equation for U (t) and solving the resulting Dyson 
series one can formally derive the Gell - Mann and Low magic formula [11] for the scattering operator 

S = T{exp(i f dt [H{t) - H })}, H(t) = e - iH ° t/h H e iH ° t/h (3.3) 

where the time ordering symbol T asks to order the latest operator to the left. Unfortunately, all of this 
is mathematically ill - defined: A tiny subsetH of the Haag - Kastler (or Wightman) axioms is sufficient to 
establish that the operator V(t) is the identity operator (up to a phase). This is Haag's famous theorem 
|llj . In other words, either there is no interaction or the magic formula is wrong. Indeed, (|3.3I) is ill - 
defined in perturbation theory and needs renormalisation. In order to avoid the implication of Haag's the- 
orem one can, as a regularisation, break translation invariance of the Hamiltonian in an intermediate step 
by multiplying its density by a function of compact support and then extend the support to infinity. This 
is also the technique underlying causal renormalisation theory [34j . 

7 The four momentum squared operator should have a pure point point spectrum which is separated from the continuum. 

8 More in detail, in order to derive the LSZ formulae one needs the LSZ asymptotic conditions which state that the matrix 
elements of the interacting field between vector states in TL+ and TL- respectively converge to those of the free field. If asymptotic 
completeness holds, this is just weak convergence on TL which is implied by the strong convergence of the Haag - Ruelle theory. 

9 It is sufficient to retain the 1. uniqueness of the vacuum and 2. spatial translation invariance of the Hamiltonian (part of 
the Poincare algebra) and 3. spatial translation invariance of the vacuum. 
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It transpires that in ordinary QFT the scattering matrix is directly related to the time ordered n - point 
functions. As we are interested in applications to quantum gravity, we are in a somewhat different situation 
because we do not have the axiomatic framework of ordinary QFT at our disposal which relies on the metric 
considered as a background field. However, one can consider a Born - Oppenheimer type of approach with 
a representation of 21 in which the three - metric operator q acts by multiplication (see [35] for first steps 
towards a technical implementation). Then, at fixed metric argument of the vector state under consider- 
ation one can consider the resulting matter part of the Hamiltonian and apply the techniques of QFT on 
curved (in this case ultra - static) backgrounds [7] and the corresponding perturbation theory [36] in order 
to define scattering theory for matter. In particular, LSZ type of formulae then again apply. To define 
scattering theory for gravity in a background independent way one should consider background independent 
semiclassical states which are concentrated on a given three geometry and extrinsic curvature and identify 
their excitations with scattering states, see |35j . 



3.2 Path Integral for n — point functions 

We are thus interested in the time ordered n— point functions. More in detail, suppose we have a repre- 
sentation of the * algebra 21 generated by the elementary fields Q A , Pa (or the corresponding C*-algebra of 
Weyl elements) on a Hilbert space 7i which supports the Hamiltonian H of the (conservative) system. We 
will assume that H is bounded from below and has at least one normalisable vacuum i.e. a unit vector 
state of minimal energy E = inf (a(H))) which is a cyclic vector for 21. Without loss of generality we redefine 
H such that E = 0. Consider the Heisenberg picture operators Q A (t) = e~ lHt / h Q A e lHt / h . As motivated 
in the previous subsection, we are interested in the time ordered n— point functions 

T A '- A "(t u ..,t n ) :=<n,T{Q A '(t l )..Q A -(t n )}n> (3.4) 

For n > 1 and pairwise distinct times we have 

re-l 

T A ^(h,..,t n ) = J2 II -**(*+!))] WMl) " Mn) (Ul),~Mn)) ( 3 - 5 ) 

ttGSu k=l 

where we have defined the unordered Wightman functions 

W Al - A »(h,..,t n ) :=<n,Q M (t 1 )..Q A "(t n )n> (3.6) 

We should allow for more general operator insertions but Q A (t) contains information about Q A (0) = 
[H,Q A ]/(ih) which knows about Pa, hence any scalar product between vector states in the dense sub- 
space 210 can be approximated by linear combinations of the functions (|3.6h . Conversely, given suitable 
positivity requirements on the Wightman functions and their transformation properties under time trans- 
lations we can reconstruct TC, f2, H via the GNS construction. The latter arises via Stone's theorem from 
the fact that we can define a strongly continuous unitary group of time translations. 
Using HQ, = we may write 

W Al - An (h,..,t n ) :=< n,e i( - t +- t ^ H ^Q Al e iH{tl - t2)/h Q A2 ..Q An e iH{tn - t - )/h n > (3.7) 

for any t±. By inserting resolutions of unity it follows that for suitable choices for ipi,ipf and times ti,tf we 
are interested in the matrix elements 

< V/, U(t f - ti)iPi > H , U(t) = exp(itH/h) (3.8) 

of the evolution operator between initial and final vectors prepared at initial and final times ti, tf respectively. 
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The path integral substitute for (|3.8p is heuristically obtained by skeletonisation of the time interval \U, tf] 
followed by insertions of unity in terms of generalised position and momentum eigenvectors respectivelM 10 !. 
Specifically, assuming that 7i is a representation in which the operators Q A act by multiplication, for time 
steps e = (tf — ti)/N and integration variables Q n := Q{U + ne), P n := Q(ti + ne) we obtain formally 

. N N N 

<lPf,U(t f -ti),^>= {\{[dQn]}{\{[dPn]}JJ{Qn)^i(Qo)[\{ <Qn,e UB ^P n > < Pn, Qn-1 >] 

J n=0 n=l n=l 

n (3 ' 9) 

where formalhl 11 ! 

[dQ] := Y[dQ A , [dP] := \{dP A (3.10) 

A A 

The assumption is now that as N — > oo we may approximate 

< Q„, e ieff / R P n >«< Q n , P n > e *H(Qn,Pn)/h (311) 

which can be heuristically justified by expanding the exponential in powers of e, ordering momentum and 
configuration operators to right and left respectively and neglecting all higher h corrections. For certain 
Hamiltonian operators of Schrodinger type one can actually prove (|3.1ip (Trotter Product formula [1]) but 
in general this is a difficult subject. Making this assumption and using the position representation of the 
momentum eigenfunction 

n expHQ^) (312) 
A ^ 

we obtain formally 

„ N N 

<^ f ,U(tf-U),A>= / iT[[dQn]} {Y[[d(P n /Vto)]} t^QnjMQo) X 

J n=0 n=l 

e N O a -O a 

exp(- i | ^ n ^ n - l P An ] - H(Q n ,P n )}) (3.13) 

n=l A 

One now takes N — > oo and formally obtains 

<rl> f ,U(t f -U),il>i>= J [DQ] [DP/sqrt2^ f {Q(t f )) MQ(U)) exp(-^£ / dt {[J] Q A P A ] - H(Q, P)}) 

(3.14) 

where 

[DQ}= II Y[dQ A (t) (3.15) 

te[u,t f ] A 

and similar for [-D-P]. If the Hamiltonian is at most quadratic in P then one can formally perform the 
momentum integral. As an example, consider a Hamiltonian of the form 

H(Q, P) = ho AB (Q)P A P B + V(Q) (3.16) 



10 This assumes that the operators Q, P obey the canonical commutation relations. For more general algebras generalised 
eigenvectors may not exist because e.g. momenta do not commute with each other. In this case one must use different resolutions 
of the identity. We will here assume that 21 obeys the CCR, CAR and more general algebras can be treated analogously. 

n There is no Lebesgue measure in infinite dimensions. However, if the Hilbert space 7i is rigorously defined as an L2 
space with a probability measure on a distributional extension of the classical configuration space, then (|3.10[) can be given a 
meaning. We will not consider these issues for our heuristic purposes and confine ourselves to drawing attention to the missing 
steps involved. 
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Examples of such Hamiltonians are for example the Hamiltonian constraint in General Relativity (neglecting 
the issue of gauge invariance for the moment) where the non trivial "supermetric" G AB {Q) is the Wheeler 
- DeWitt metric and the potential V(Q) is related to the Ricci scalar of the three metric Q. (In)famously, 
neiher G nor V are positive definite so that the Hamiltonian is not bounded from below in General Relativity. 
In any case, for Hamiltonians of type (|3.16p we can formally perform the Gaussian integral and obtain 



< 



rP f ,U(tf-ti),i>i >= W [DQ] Mdet(G)|] ^/(Q(i/)) MQ(U)) exp(-i /*' & {[^(G~ 1 )abQ A Q B -V(Q)}) 

(3.17) 



where M is an (infinite) numerical constant (a power of 2tt and h) and 



h/dit(G)]= J] ^ditM (3.18) 

te[ti,t f ] 



is the functional determinant of the supermetriq 12 !. 

Notably, if G is a non trivial function of Q then it is not true that 



<^ f ,U(t f -t i ),ij i >=M I [DQ] ip f (Q(t f )) ipi(Q(ti)) exp(--S[Q,Q;[ti,t f ])) (3.19) 



with the classical action 



S[Q,Q;[ti,t f ]]:= V dt L(Q,Q), L(Q,Q) = ^(G~ 1 )abQ A Q B ~ V(Q) (3.20) 

Even worse is the case that the momentum dependence of the Hamiltonian is higher than quadratic so that 
the momentum integral can no longer be performed exactly. In that case one can at best perform a saddle 
point approximation or one has to rely on perturbation theory. We see that the correct path integral in 
general is over the phase space and involves the Hamiltonian action and not only over the configuration 
space involving only the Lagrangian action, so we will stick with (I3.14|) in what follows. 

We still must provide a path integral formulation for the n— point functions. However, this is is easy by 
noting that 

n 

W Al " An (h,..,t n ) = H / [dQ k ] <n,U(t + -h)\Q 1 >Qf 1 x 
fc=i ^ 

n-l 

[JI < Qk, U (tk ~ t k +l)\Qk+i > Qt+i 1 } < Qn\U(t n - t-)\Sl > (3.21) 



k=l 



where Q A \Q >= Q A \Q > was used. Combining (|3.2ip with (|3.14[) results in (for t + > t\ > ..t n > i_) 
W Al - A "(ti,..,t n ) = J [DQ] [DP/V2^MQ(t + )) fi(Q ( t_)) x 

/*+ n 
dt{[Y^Q A PA]-H(Q,P)}) \{Q A Ht k ) (3.22) 
- 1 k=i 



where 



[dq]= n ( 3 - 23 ) 

te[t-,t+] A 

and similar for [DP]. 



12 In fact there is a sign factor involved which accounts for the signature of G. Equation (|3.17[) is only correct if the signature 
of G does not depend on Q. 
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It is worth mentioning that in a rigorous setting [TJ [H] one does not really consider matrix elements of 
the unitary operator U(t) = exp{itH/K). Namely, consider the analytic continuation t k \— * itj~ for t k > 0, 
that is, the Schwinger functions 



S Al - An (h,..,t n ) :=W Ai - An {itx,..,it n ) (3.24) 

These are correlators of the e tH / h Q A e - tH / h anc l now the same formal manipulations as before lead us to 
consider the contraction semi - group t i— ► V(t) = exp(—tH/h), t > 0. One now obtains instead of (13,221) 
the formula 

S Al - A -(t u ..,t n ) = J [DQ] [DP/V2^MQ(t + )) fi(Q ( t_)) x 

i /•*+ n 
exp(--/ dt{[i^2Q A P A ] + H(Q,P)}) HQ Ak (t k ) (3.25) 

*J t— A 7 .1 



fc=l 



For Hamiltonians of the form (|3.16j) with positive definite G, V (subtract the energy gap if necessary) the 
formal Gaussian integration now gives 

r i n 

S Al " A »(t 1 ,..,t n )=M / [DQ] [y/deb(Gj] n(Q(t+)) n(Q ( t_)) exp(--S E [Q,Q; [*-,*+]]) 11^^) 
J k=i 

(3.26) 

with the "Euclidian" action 

S[Q,Q;[ti,t f ]]:= V dt L E (Q,Q), L E (Q,Q) = \{G- 1 ) AB Q A Q B + V(Q) (3.27) 

The path integral (13. 26ft has better chances to be rigorously defined because the "measure" has a damping 
factor rather than an oscillating one and so in the rigorous setting one defines (|3.22|) by backwards analytic 
continuation of (|3.26p (when possible )EE Equation (|3.26p (when it can be proved) is called the Feynman - 
Kac formula [DEIS]. in what follows we therefore consider the Euclidian point of view. 

In order to avoid any infinite constants we divide the contraction matrix by 1 =< f2, £1 >=< 0, V(t+ — 
t-)£l > and obtain formally 

sM t = J[DQ] [DPMQ(t+)) exp(-i//_ + dt{[iY. A Q A PA]+H{Q,P)})Y[ n k=l Q A K{t k ) 

n[h, - ,tn) j [DQ] [DP]Q(Q(t + )) n(Q ( t_)) ex P (-i//_ + dt {{i Q A Pa] + H(Q, P)}) 

(3.28) 

Even if one cannot inegrate out the momenta in general, formula (|3.28p reveals that what we are interested 
in is the measure formally given by 

1 



dfJt(Q) := - [DQ] exp{-S E [Q]/h) n(Q(t+)) (3.29) 

Zj 



where 

i /■*+ 



exp(-S E [Q]/h) := J [DP] exp(-i^ + dt {[% ^ Q A P A ] + H(Q, P)}) (3.30) 
is the exponential of the effective Euclidian action and 

Z := j [DQ] exp(-S E [Q]/h) Q(Q(t+))n(Q(tJ)) (3.31) 



13 It is worth mentioning that in teh axiomatic framework of local quantum physics [11] on Minkowski space the Schwinger 
functions are automatically symmetric although the Wightman functions are not which is a consequence of the locality axiom 
(bosonic operator valued valued distributions supported at spacelike separated points commute) and analyticity. In GR one 
does not expect to construct a Wightman QFT due to background independence which is why we insist on > t^+i- 
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is the partition function. None of the three quantities [-DQ], Se, Z exists but in fortunate cases their 
combination can be rigorously defined as a measure on a suitable distributional extension of the space of 
configuration variables Q. The measure \i is known if we know all its moments or equivalently its generating 
functional 

X [j] := J dfi(Q) e lE ^- dt ^ QA ^ (3.32) 

from which the moments follow by (functional) derivation at zero current j. 

The apparent drawback of the these formulae is that they involve the exact ground state O of the 
interacting Hamiltonian H which is difficult if not impossible to compute analytically. However, and here is 
where the Euclidian formulation again is helpful, notice that so far the choices for t± were arbitrary except 
that tk £ [t~, t+], in particular, in the original correlator the dependence on t± is through e~ t+H £l = and 
e t ~ H Q = £1. Now suppose in addition that exp(-tH) for t > has a positive integral kernel, i.e. maps 
a.e. positive functions to strictly positive functions which is usually the case. Then it follows from [37] that 
E = is a simple eigenvalue and the unique (up to a phase) ground state is a strictly positive function. 
It can be obtained from any a.e. positive Oo £ TC via the strong limit 

O := lim 6 "° (3.33) 

It follows that by taking the limit t± — > =too we can replace fi by Qq in f|3.29[) - (|3.32[) because the factors 
of ||e _i ^Oo|| cancel in numerator and denominator. We will assume this to have done for what follows. 
Remarkably, the choice of the reference vector Qq is rather arbitrary. 



Having justified the replacement of by f^o i n the Euclian regime, we analytically continue the time 
parameter backwards to define the time ordered n-point functions and thus the exponential becomes a pure 
phase. 



3.3 Unfolding the Reduced Phase Space Path Integral 

We would like to rewrite the path integral over the reduced phase space coordinatised by the chosen true 
degrees of freedom in terms of the unconstrained phase space. This is of course standard, see e.g. [6], but 
we review this procedure here for the sake of completeness. It is, however, a rather involved procedure. 



3.3.1 Preliminary Results 

The virtue of the gauge fixing conditions G is that the system C := {S, F, G} is now a total second class 
system so that one can treat all constraints on equal footing. We will do this first in the adapted system of 
Darboux coordinates (Q , Pa), O^jT/i) {x a ,y a ) which is related to the original system (q a ,p a ) by a (local) 
canonical transformation and then show that the resulting expression is actually invariant under canonical 
transformations. 

Theorem 3.1. 

Let C = {Ca\ be a second class system of constraints on a phase space with canonical coordinates z 1 and 
symplectic structure uj on the unconstrained phase space V. Denote the constraint surface by T := {m £ 
T; Cj^{m) = V ^4} which is a submanifold ofT. Consider an embedding J : T — > T with J(T) = F where 
F is a model manifold of with coordinates x l for F . 
i. 

uj := J*lo is a symplectic structure on F. 
ii. 

Let £1* be the degenerate symplectic structure on F defined by the Dirac bracket corresponding to C . Let 
f,g £ C 1 (r). Then J*({f,g}*) = {J*f, J*g} A where {., .} A is the Poisson bracket associated with uj. 
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Hi. 

The relation between the Liouville measures [II and (ll onT and T respectively is 

Ax [J*/] = /U L [Vdet({C,C}) 5(C) f] (3.34) 

for any measurable function f . 

We note that the right hand side of (|3.34p does not make any reference to the chosen embedding J. 

Proof. 
i. 

Obviously du = J*dVt = establishes closure. Non - degeneracy follows from the fact that J has maximal 

rank. 

ii. 

Let 

Mab := {C A ,C B } (3.35) 

thenB 

[oj*) ij = L0 IJ + (M- 1 ) AB u) IK oj JL C AjK dC BjL (3.36) 
where u IJ u jk = <$x- Using that {/, g} = lo ij f t jgj and 

{J*fi J*g} A = & j (J*f),i (J*g),j = & 4 f,i Jfj g,j (3.37) 

with u l3 ijjjk = Si, we see that the claim is equivalent to 

LU^Jpfj = [lo*] ij (3.38) 

on M. To verify (13. 38ft we notice that a A := (C A j), <7j := (ujjjjf) is a linearly independent set of one 
forms on M and it suffices to check (|3,38p in this basis. From J*C A = for all A we immediately have 

JpA,I = (3.39) 

on M and by construction of the Dirac bracket it is not difficult to see that contraction of (|3,38p with a a 
results in zero on both sides. Contraction with OiOj results in the identity 

& 1 J,k Jj a U a jJ = U M &lj = tiji 

= [io*] IJ a iT a 3 j = oj ij an ojj = u IJ co IK uj jl jf J L j uj lk jf jj (3.40) 
where we used (|3.39|) and 

Uij = (J*u)ij = ujlj 4 Jj (3.41) 

iii. 

Recall that for finite (2n-) dimensional systems the Liouville measure is simply [i^ := f\ n oj = Vi(u))[dz\ 
where Pf (uj) = y / det(w)) denotes the Pfaffian of the matrix ujjj. We adopt here the same formula for infinite 
dimensions, ignoring as ususal that the Lebesgue measure [dz] does not exist. Using (I3.35|) we solve the 5 
distribution in (|3.34p in terms of the embedding J which we write in the form z = (x, y) = J(x) = (x, Y(x)). 
Here x, y are separate sets of canonical pairs so that ujjj becomes block diagonal and the block matrices 
^AB, ^ij are constant. We obtain 

» L Wdet({C,C}) 6(C) f} = J [dz] ^det(u)det(M)(z)5(C(z))f(z) = f [dx] ( ^M^p /)(J(x)) 

(3.42) 

14 Our conventions are as follows: i Xs ^> + df := defines the Hamiltonian vector field \f associated to / while {/, g} := 
~Xf[g] ~ ~ixjdg = ix/*x B tJ defines the Poisson bracket. The corresponding matrix is denoted by u IJ := {z 1 ,z J }. 
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where c AB := C A>B . Here we used C A (x,y) = C A (x,Y(x)) + c AB [y — Y{x)] B + .. = c AB [y ~ Y(x)] B + ... We 
have 

M AB = lo ij C aj C Bi j = lo cv C A:C C b , v + LO^C Ati C Bij (3.43) 
Equation (|3.39p takes the form 

C Al + C AB Yf = (3.44) 



so that (|3,43p can be written 

Let us introduce the abbreviations 

then 

Consider now the matrices 
The key identity is now 



M AB = c ac c bv [uj cv + u^YfYf ] (3.45) 

Y A :=Y A , Y\:=u AB ^Y B (3.46) 

M AB = c ac c bv lo £V [5 C £ - YfY*} (3.47) 

K A := Y A Y^ k{ := Y A Y A (3.48) 



det(l — K) = det(l - k) (3.49) 
To prove this we use the identity (supposing that k has rank m) 

m 

det (l - k ) = 1 + ("!)' <f "4? k n- k n ( 3 - 5 °) 
i=i 

The same formula holds for det(l — K) just that K may have a different rank n and that summation indices 
are A rather than i. Now each term in the sum of (13.500 is a polynomial in the the traces tr(/c r ), r > with 
a coefficient that does not depend on m. However, tr(/c r ) = tr(K r ) for any r. So the only possible difference 
in the two quantities is the range of I. However, notice that 

5f fcf ..0 = Yf.Yf Y A \.Y Al = Y}\ .Y k A , Y Al ..Y Al (3.51) 

3\ 3i n h Ai Ai ii ii [Ai Ai] ii ii v j 

is completely skew in either set of indices, hence the sum anyway extends to min(m, n) only. 
We conclude with det(cj) = det((ujij)) det((u> AB )) that 



d6t( d ^ 2 t ) (M) = det((^)) det($ - #)) = det(K- - UM YfY* 



det((w ifc )) det((^ - tf)) = debfaij - uj AB YfY^) = det(c^) (3.52) 

□ 



Corollary 3.1. 

The measure /Ug on T defined by ft3.34\ ) in terms of a gauge fixing condition G, as a linear functional is in 
fact independent of the gauge fixing condition when restricted to gauge invariant functions f . 

Proof. By definition of a gauge fixing condition G for a first class constraint set {F}, it defines a section of 
the first class constraint surface (i.e. it defines a hypersurface that intersects each gauge orbit in precisely 
one point) and it can be reached from any point on the same gauge orbit. Hence any two gauge fixings G, G' 
are related by a gauge transformation ip which can be written as a composition of canonical transformations 
of the form exp(/3^'{i ? /J , .}) for real valued (phase space independent) parameters. By the first class property, 
there exist matrices L, M, N such that ip ■ = L^F U and ip ■ S% = M|? S^y + N^F^ where L, M are non - 
singulaJ^l. In matrix notation ip ■ F = L ■ F, ip-S = M- S + N-F. This can be inverted 

F = (L- 1 ) ■ ((p ■ F), S = (M- 1 ) -[(p-S)-N- (L- 1 ) • {<p • F)] (3.53) 



3 At least for /3 M close to zero. 
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By assumption, / is (weakly) gauge invariant, f{m) ~ /(</? • m) and the Liouville measure is invariant under 
canonical transformations (since the symplectic structure is), dfii(ip ■ m ) = c£//£,(m). 

We exhibit the dependence of the measure (|3.34p on G by Notice that in terms of {C} = {F, G, 5} 
we have 

/ {F»,F V } {F^,G U } {F„,Sv} \ ( {F^,G U } 

det({C,C})c;=o = det {G^,F U } {G^,G U } {G^Sjy} = det {G^} {G^C} {G M ,5 S '} 

\ {5s, F,} {5s, G^} {5s, 5s'} / c=0 \ {5s, G^} {5s, 5s'} , c=0 

= {[det({F, G})] 2 det({5, 5})} c =o (3.54) 
Using the automorphism property of canonical transformations [tp ■ f] (m) = f{ip- m) etc. and f)3.53|) we have 

/Vg(/) = / dfi L {m) S(S(m)) 5([<p ■ G](m)) 5(F(m)) \ det({F, <p ■ G}(m))\ ^/det({S, S}(m)) /(to) 

dfi L {m) 8((M~ l [(p ■ S — NL~ l (p ■ F]){m)) S(G(ip ■ m)) 5((L _1 [^ • F])(m)) | det({L -1 ^ • F, p • G}(m))\ x 

M 



Vdet^M- 1 ^ • S - JVL- V ' F], M-^tp ■ S - NlFhp ■ F]}(m)) f{tp • to) 

dfi L (m) | det(M)(m)| |(det(.L)(m))|tf(fe> • 5](to)) <J(G(<p • to)) *(fo> • F](m)) | det^IT 1 ^ ■ F V ' G})(m))\ x 



Y / det(M- 1 {< < 5 ■ 5, y> ■ 5XA/- 1 ) 2, - NL'^ip ■ F,ip- SjiM- 1 ) 7, 

-M-ifa ■ S, <p ■ F]}(NL- 1 ) T + NlF^Jp -F,ip- ^]}(^'L- 1 ) t )(to)) f(<p ■ to) 

f dfJ, L (m)\det(M)(m)\ \(det(L)(m))\S([(p ■ S](mj) 5(G(<p ■ m)) S([ip ■ F](m)) | det((X — ■ F,tp ■ G})(m))\ x 
y/det(M-i{<p -S,<p- S}(M-f)(m))/(^ • to) 

/ dfj, L (m) 6([<p ■ S](m)) 6(G(<p ■ m)) £([<p • F](m)) | det(({<^ • F, <p • <7})(m))| x 

y/det{{ip-S,ip-S}(m))f((p • to) 

■ to) <5(5(<^ • m)) 5(G(y> ■ to)) S(F(tp ■ to)) | det(({F, G})(tp ■ m))\ x 



,1/ 



y/(det({S,S})(<p-m))f(<p ■ to) 
= M/) ( 3 - 55 ) 

where in the third step we used that Poisson brackets with L, M, N do not contribute since the 5— distributions 
have support at (p ■ F = <p ■ 5 = 0, in the fourth we used the first class property and again the support of the 
5 distributions, in the fifth we cancelled the determinants of the matrices L, M, in the sixth we exploited 
the Poisson automorphism property of ip as well as the invariance of the Liouville measure and in the last 
we performed a trivial relabelling. □ 

The statements of theorem 13. II and corollary 13.11 show that the measure \xq (|3.34p is the correct extension 
to the full phase space of the pull - back measure defined by a gauge fixing condition and that correlators 

among gauge invariant functions are actually independent of the gauge fixing condition. For instance, in 

(G) 

terms of the gauge invariant observables Oj , where we have exhibited the dependence on G, we have 

HG,[0f ] ] = ^G[0f ] ] for any G' = (p ■ G. 

This can also be understood geometrically: Given two gauge fixing conditions G, G' we obtain log = 
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JqUJ, ll>g> = Jq/UJ from the corresponding embeddings Jq : M —> Mq, Jqi : M — > Mqi. Now clearl\f^l 



M V . G 



{to e M\ S(m) = F(m) = cp*G(m) = 0} 

{to, G M; M^m^S^M) - N(m)L- 1 (m)F(ip(m))]L- 1 F(<p(m)) = G(<p(m)) = 0} 
{to G M; S(tp(m)) = F((p(m)) = G{ip{m)) = 0} 
{ip-\ip(m)) E M; S(<p(m)) = F(<p(m)) = G(<p(m)) = 0} 



(3.56) 




(3.57) 



(3.58) 



Remark: 



The fact that ip is canonically generated by first class constraints featured crucially into this argument. This 
has the following relevance: 

Suppose we are given a system which as gauge symmetry has spatial diffeomorphism invariance in D spatial 
directions. Suppose that the field content consists, possibly among other things, of GR minimally coupled 
to D scalar fields <j>i, ..,(f> n . From the curvature of the metric and higher derivatives we can also form D 
algebraically independent scalars Ri,..,R n . Suppose that at least locally they define a coordinate system 
so that x i— > tfi(x) := (<f>i(x), .., (j> n {x)) and x i— ► R(x) defines a (local) diffeomorphism. Pick any fixed 
diffeomorphism ip®. Then both G = <j> — (fo and G' = R — <po are bona fide gauge fixing conditions. However, 
there does not exist any canonically generated diffeomorphism ^ = exp({J d D x ^ a {x)C a {x), .}) with phase 
space indendent £ such that <p% • G = G' . The reason is that the spatial diffeomorphism constraint does 
not mix field species. It is true that we can find a phase space dependent function £[(/>, R] defined by 
(pg o <j) = R such that • G]^—^ = G' , however, due to the phase space dependence of £ it is not true that 

[v^°]^ = | = P^ - The latter is also a canonical transformation with generator f d D x £ a C a but it does not 
generate the searched for field dependent diffeomorphism, provided it exists at all. Notice that corollary 
13.11 remains true for field dependent £, just the matrices L, M, N look different, this is not the point, the 
point is that it is not clear that a canonical transformation exists which maps <j> to R. It may therefore be 
true that gauge fixings seperate into equivalence classes depending on whether such phase space dependent 
gauge transformations exist or not. If that was the case, then it would not be true that the measure ()3.34l) 
as a linear functional on gauge invariant functions would be independent of the gauge fixing condition, it 
would depend at least on the equivalence class. 

3.3.2 From reduced Darboux coordinates to unreduced Darboux coordinates 

In order to combine the results of sections 13.21 and 13.3.11 we notice that the parameter manifold M (which 
is the same for any gauge fixing) can be identified with the manifold equipped with Darboux coordinates 
{Q a ,Pa}- These are adapted to our choice of G such that F = G = or equivalently F' = G = can 
be solved for {t^Tfy} in terms of {(2 a ,1U} which also defines the embedding Jq- In particular, if / only 
depends on {Q , Pa} then we can form our preferred observables O ~ and due to the identity / = JqO ~ 



16 In abuse of notation we write ip ■ m = <p(m) i.e. we identify the action of the exponential map with the corresponding 
diffeomorphism. 



we find from f|3. 34[) 



Al[/] = Mg[0? ] = fl><p-G[Of] 



(3.59) 
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where corollary 13.11 was used. Of course, for practical calculations the precise expression for O ~ in terms 

of Q A , Pa, ft/I is rather cumbersome to use. However, due to the 5 distribution 5(G) involved in fiQ 
obviously 

(i L [f}=fi G [0f]=f J/G [f) (3.60) 

so that we can drop the gauge invariant extension under the path integral at the price of having to to keep 

the G dependence in fj-cG because / is not gauge inavariant so that corollary 13.11 does not apply. Even if 

(G) 

we keep O ~ rather than /, still the G dependence does not disappear because while we can drop it from 

(G) 

fiG, it remains in ~ which is a specific type of Dirac observable which uses the structure G. This is in 
accordance with what we said in the introduction. 

We are now ready to extend the reduced Darboux coordinate phase space path integral of section 
13.21 to all Darboux coordinates: The Liouville measure used there is precisely given by jli because in 
Darboux coordinates det(o>) = 1. Furhermore, for our choice of gauge fixing G^ = + r^(t) and 
FL = tt^ + h'^^jQ, P) we have | det({F' , G})\ = 1 and since S' a = zy, = (x cr ,y CT ) in Darboux coordiantes 
are canonical pairs we have det({S",5"}) = 1. It is therefore trivial to write the generating functional of 
n— point functions as a path integral over the entire phase space by simply using formula (13.341) at each 
point of time 

Z\j) := J [DQ DP D<t> Dtt Dx Dy] 5[G] 5[S'] S[F'] | det[{F',G}| y/det[{S' , S'}] n (Q(+oo)) ft (Q(-oo)) 
exp(-4 / dt {[Y,Q A Pa] ~ H rcd (Q,P;t)}) e*/« dt jA (t)Q^t) (3 _ 



where for instance 



and 



6[F'} = [] 6(F'(t)), det[{S', S'}] := ]J det({5'(t), S'(t)}) (3.62) 



S A - A "(t u ..,t n ) = r n [ ] J=0 (3.63) 

Here we have explicitly kept det({5", S'}) = 1 because we will see that (|3.61|) is covariant under changing to 
equivalent constraints. To remind the reader, we recall that the possibly explicitly time dependent reduced 
Hamiltonian is given by 

H rcd (Q, P; t) = rUfyti^ = r(t), Q, P)) (3.64) 

where F^ at S = or equivalently S' = was brought into the equivalent form F' = ir^ + h'((f>, Q, P) which 
motivated the use of a gauge fixing of the form G^ = T^(t) — ^ . 

Formula (|3.6ip achieves the goal to extend the reduced phase space path integral to the full phase space, 
albeit in the specific, local Darboux coordinates that were picked by motivations from quantum theorM 17 ! 
and the constraint structure of the theory and in terms of the convenient equivalent constraints S',F'. 



lr Due to the second class constraints, the use of such coordinates is mandatory because otherwise the representation theory 
of the reduced symplectic structure becomes too difficult. 
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3.3.3 Restoring the Original Canonical Coordinates and Constraints 



The next step will be to restore the original Darboux coordinates (q a ,p a ) as well as the original constraints 
S, F rather than S',D' . To that end we notice the identity 

dt [P A Q A -H rcd (r;Q,P)} 

dt \P A Q A -f^{T-Q,P)] 

dt [P A Q A + tt^ + Tr^r" - 01 - + h^r- Q, P)}] 

dt [P A Q A + tt^ 1 + 7r^G M - T»Fp\ (3.65) 

Since the path integral is supported at = — ^ = 0, F'^ = ir^ + h' = 0, := zy, = (x a , y a ) = we 
can rewrite ()3.65j) under the integral in the form 

J dt[P A Q A + n^ + y a y a ) (3.66) 

Now for certain, phase dependent, non singular matrices M, N we have F' = M ■ F, S' = N ■ S. But then 



S(F')\dei({F',G})\=S(F)\det{{F,G})\, 6(S')y/det({S',S'}) = 6(S)y/det({S,S}) (3.67) 

is covariant under change to equivalent constraints. This allows us to immediately restore the original 
constraints in (|3.61j) altough everything is still written in terms of the unreduced and adapted darboux 
coordinates (Q, P), (<j), ir), (x, y). However, that system of coordinates originates from the original system of 
canonical pairs (q a ,p a ) by a a canonical transformation [6]. Accordingly, by applying the inverse canonical 
transformation a, we can restore the system of coordinates (q,p) which leaves the Liouville measure in (|3.6ip 
invariant, which leaves (|3.66j) invariant up to a total differential which we assume to vanish at t±, which 
reexpresses F, G, S in terms Darboux coordinates in terms of the original coordinates and finally is covariant 
with respect to the Poisson brackets involved because e.g. a({F, G}) = {a(F), a(G)}. Therefore, (|3.61l can 
be rewritten as 

m ■= f [DqDp] S[G] S[S] 6[F] \det[{F,G}\ ^det[{S,S}] n (Q[q,p\(+oo)) n Q (Q[q,p](-oo)) x 

/ dt [VVpa]) e*/* dtjA(t)Q A [q, P ](t) (3. 68 ) 



ex Pl -- 



Notice that the effect of the reduced Darboux coordinates did not completely disappear: The initial and 
final state depend on Q as well as the exponential involving the current j. But Q = Q[q,p] may be a 
complicated function of the original canonical coordinates q,p. 



Remark: 

Notice that at this stage we can formally get rid of the gauge fixing condition in (|3,68p by the "Fadeev - 

Popov trick" if we pay a price: As we have already remarked before, due to the presence of 5[G] we may 

(G) 

replace everywhere the non gauge invariant Q by Oq . Then the exponent, as a symplectic potential and 
the measure [Dq] [Dp], which is formally the Liouville measure associated with u, is gauge invariant since 
gauge transformations are canonical transformations. Also {S, S}, F are weakly gauge invariant due to the 
first class property and since canonical transformations preserve Poisson brackets. Hence, after the gauge 
invariant extension of Q, the only non - gauge invariant ingredient of the integrand of Z[j] in (I3.61|) is 
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S[G]. In fact, ap(G) = G — (3 where we have introduced the gauge transformations ag := exp(/3 M {i^, .} 
which, since the constraints F' are Abelian, have the Abelian group structure ap o a pi = ap+p>. Since 
the remaining ingredients are all gauge invariant, we may replace G by ap(G) for any (3. Now extend both 
numerator Z[j] and denominator Z[0] by the infinite "gauge volume" constant f[D/3] with "Haar measure" 
[D0\. We can then trivially integrate out the 5[ap(G)] = 5[G — 0] and find 

X[j] ■= ^ 
xm Z[0] 

Z[j] := J [Dq] [Dp] S[S] S[F] y/det[{S,S}] ^ (og ?) [g,p](+oo)) Q (O^ \q,p](-oo)) X 



- f d tJ2Q a P a )e^ dtjAit) °>' Pm (3.69) 



The price that we have to pay is that we have to replace Q A by Oq5 =: Oq a e ^ c - which is an even more 
complicated function of q,p. This makes this method of getting rid of the gauge fixing condition useless in 
practice. The only exception is when we consider zero current j = and restrict to gauge transformations 
that are the identity in the infinite past and future. Then the gauge invariant extension in the argument of 
Qo and more general boundary states is not necessary and the argument just displayed goes through. The 
restriction to such asymptotically trivial gauge transformations however means that we set the constraints 
to zero by hand on the kinematical Hilbert space. 

3.3.4 Restoring the Lagrangian 

The idea is now to exponentiate the constraints and to rewrite the total resulting exponent in terms of 
the classical action. Thus we introduce Lagrange multipliers for all first class constraints and for all 
second class constraints and write 

xm ■■= m 

Z[0] 

Z[j] := / [Dq] [Dp] [DX] [Dfj] 6[G] |det[{F,G}| vMettfS.S}] n (Q[g,p](+°o)) n (Q[q,p](- oo)) x 



exp(-l f dt [Y,q a Pa -}Z^F, ~ 5> E &0]) dt ^ QA ^ (3.70) 

J R a fj, E 

The final task is to remove the secondary constraints so that the action appears in covariant form after 
Legendre transformation, i.e. with primary constraints only. The technique for doing this is well known 
|12j . We will not treat the general case with secondary constraints of arbitrary high but finite order N (i.e. 
one has secondary, tertiary, .. N-ary constraints). For a systematic classification of such systems and a 
convenient choice of basis of those constraints see e.g. [38] and references therein. Here we pick a typical 
situation of particular interest for General Relativity. The general case will be even more complicated with 
even more modifications to the measure than we encounter below: 
We assume that the canonical Hamiltonian has the following structure 

H = H' +v f ■ F^ +v s - S^, HV) = -qfPW (3.71) 

As the notation suggests, F^ is a vector with entries consisting of primary first class constraints while 
is a vector with entries consisting of primary second class constraints. The vector F^ is related to a set 
of secondary second class constraints that we will derive shortly. Usually F^ , simply express the fact 
that the momenta pj,p s conjugate to qj, q s respectly vanish because the Lagrangian does not depend on 
the velocities Vf,v s . It is also usually the case that F^> does not depend on q/,Pf,p s but on q s . Thus we 
assume that (we do not denote indices) 

{F (i) )jF (i) } = (iXi), <?(!)} = = {F^,F^} = (3.72) 
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while the terms not proportional to , of 

S (2) . = {i H2) ;(S (l) } {3 73) 

defines a vector of secondary second class constraints. We abuse notation by identifying that vector with 
(|3.73|) in order not to have to introduce indices. Notice that 

{H, F«} = {H' , F^} = F® (3.74) 

thus justifying the name secondary first class constraint. 

It is also often the case that the F^ 2 ) close on themselves, that is 

{F( 2 ),iH2) }oCj p(2) ( 3. 75) 

which we also will assume. Correspondingly, 

{H,F^} = {H' ,F®} + {v s ■ SW,FW} oc F( 2 ),S( 2 ) (3.76) 
does not produce tertiary constraints. These assumptions imply by the Jacobi identity that 

{F (i) i5 (2) } = -{P®,{S( 1 \fU}}-{S( 1 \{f( 1 \fW}} = (3.77) 

Finally 

{H,S^} = g f ■ {fM,sM} + v s ■ {S (1 \S^} (3.78) 

and we assume that the matrix {5''- 1 \ S^} is invertible. Hence the Dirac algorithmus does not produce any 
tertiary constraints and the velocity v s must be fixed in order to equate (|3.78p to zero. Accordingly the 
Hamiltonian becomes 

H = qf • [F& - {pW, 5( 2 >}[{5W, SW}}- 1 ■ S«] + v f ■ F« =: q f ■ F^ + v f ■ F« ( 3 .79) 

a linear combination of first class constraints. Thus in terms of the previous notation, the first class 
constraints Fn comprise F^ , F^ , the second class constraints comprise SM,S® and finally the primary 
constraints C, comprise 

This is a simple but non trivial situation often encountered in concrete models and this concrete form 
now enables us to explicitly carry out the steps outlined in [12] : In (|3.65p by an obvious change of notation 
we write 

Z[j] := [ [Dq] [Dp] [DXi] [DX 2 ] [D^] [D^] 5[G] |det[{F,G}| ^Jdet[{S,S}] fi (Q[9,p](+oo)) O (Q[g,p](-oo)) 



exp( _I f dt {{i(52 q a Pa - Ai • F« - A 2 • F< 2 > - mi • S« - ^ • S^)]}) e*/" *Mt)Q A l<u>)® 



[Dq] [Dp] [DXi] [DX 2 ] [Dm] [D^] 5[G] |det[{F,G}| Vdet[{5,5}] O (Q[g,p](+oo)) n (Q[q,p](-oo)) x 
exp(-~ f dt {[i(Y^q a p a - Ai • F« - A 2 • F^ _ ^ . 5(1) - m ■ S (2) )]}) e*^ M ^t)Q A [i,Pl(t) ( 3 . 80 ) 

where in the second step we have shifted the integration variable m in order to absorb the contribution 
M-SW=FW-p). 

We now perform a canonical transformation with generator /i 2 ■ S^/X' 2 at each time t £ which 
we assume to become the identity at t±. Here X' 2 is the unique component of A 2 such that A 2 • {F^ 2 \ S^} = 
A2<S( 2 ) modulo terms proportional to F^\ S^. This transformation preserves the Liouville measure, the 
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symplectic potential J dtp a q a and £lo(Q[q,p](+oo)) Qo(Q[q,p](— do)) e l ^» dt iA(t)Q A (t) s [ nce [ n this example 
under consideration Q is among the coordinates independent of q s . If we assume that G = GW,G( 2 ) do 
not involve (; s then also G is invariant. Indeed, one can choose G^,G^ to be gauge fixing conditions 
on q"p := qj and qp respectively which are both independent of q s = q^ so that qs , qP , qf\ qf \ Q A 
comprises a complete system of configuration coordinates. Furthermore, clearly F^\S^ are invariant. 
But, denoting the canonical transformation by a, we have 

a(X 2 ■ F&) = A 2 • - fi 2 ■ - _L^ 2 . {S \!0} ■ ^ + 0(/i 3 2 , F« S®) 

a(/i 2 • S&) = n 2 ■ + • {S 1 , S®} ■ M2 + 0(t4) (3.81) 

Since integrating again over Ai,^i enforces F^ = = we can drop terms proportional to F^\ . 

Next, a({F, G}) = {a(F), a(G)} = {a(F), G} under the assumptions made. This will in general depend 
non trivially on \x 2 through a(F^). Likewise a({S, S}) = {a(S), a(S)} will in general depend non trivially 
on ju 2 through a(S ( - 2 ^). Consider | det({-F, G})\, det({5, S}) as expanded in powers of [i 2 - Also, since 

a(X 2 ■ F® + M2 • S^) = A 2 • F^ + ^-,u 2 • {5«, S&} • M2 + 0(/xij) (3-82) 

let us power expand a(5(F))5(S)) around that quadratic term and perform the Gaussian integral. Since 
det({S, 5}) = det{{S^ , S^}) 2 this yields 



Z[j] = J [Dq] [Dp] [DXi] [DX^ 1 / 2 ] [D/z a ] 5[G] | det[{F,G}| {/det[{S,S}] V fi (Q[?,p](+oo)) O (Q[g,p](-oo)) x 
exp(-i / & - Ai • F« - A 2 • F< 2 > - Ml • S«)]}) ^^WQ A [^1W (3.83) 

where the local factor V accounts for the additional contributions just mentioned. The integral over Pf = 
F^',p s = produces <5(Ai — Vf)8(pL\ — v s ) and cancels the integral over X\,pL\. Denoting {q' a } = 

{qf,qf\Q A }, {P' a } = Uf\qi 2 \Q A } we are left with 

Z[j] = j [Dq][Dp'] pAalA'// 2 ] 5[G] \det[{F,G}\ tydet[{S, S}] V Q (Q[q, p](+oo)) O (Q[g,p](-oo)) x 
exp(-i f dt{[i(y / q' a p' a -X 2 -F^})e i ^dtjA{t)Q A [ q ,p)(t) (3 84) 

Now, by definition (see also sectoion [2.ip . solving p' a = dL(q'^ , qj, q s ; v'@) / dv' a for v' a yields v' = u' a {q 113 , q s ;p'g) 
and 

^ = qf^Hq^q'";^) = [v' a p' a - L(q'P ,q f ,q s ;v'P)] v ,= u , (3.85) 
As is well known, the inverse of this Legendre transformation is 

L = [v' a p' a - H' {qi\q fl q s -pi p )\ pl=dL/dv , (3.86) 

Therefore a saddle point expansion about the extremum p' = dL/dv' of the exponent in (|3.84p yields 



Z[j] = J [Dq] [DX 2 \X' 2 \ 1 ^ 2 ] 8[G] |det[{F,G}| ^/det[{S,S}] V' n (Q[q](+oo)) n (Q[q](-oo)) 



x 



exp(^ / dt L{q f = X 2 ,q s ,q' b ;v' b )) e l f* dt w(*)<? A M(*) (3.87) 
h Vr 

where V'/V accounts for the additional modifications that come from the saddle point approximation. 
Notice that a possible dependence on p in Q A prohibits the saddle point approximation beyond its zeroth 
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order term! Assuming that Q is independent of p as it is the case in this example under consideration and 
assuming that really is a coordinate condition on qj and noticing that F, S, G^ 2 \ V do not depend on 
qj we have after relabelling A2 — > qj 

Z\j] = f[Dq] [W^ 2 ] 5[G®] |det[{F,G}| ^/det[{S,S}} V n (Q[q](+oc)) O (Q[g](-oo)) x 

exp(^ / dtL{q f ,q s ,q' b -v' b ))e i ^ dtjA ^ QA ^ (3.88) 

which is our final result. 

To summarise: 

The path integral can be brought into a form only involving a configuration integral and the exponent of the 
covariant action, but there is a non trivial measure factor depending on S, F, G, V' which accounts for the 
correct implementation of the dynamics. Missing that factor means quantising an entirely different system. 
The measure is not covariant with respect to the Lagrangian symmetries, however, by construction it is 
covariant with respect to the Hamiltonian symmetries generated by the first class constraints [22]. As is 
well known from classical Noether theory, these two symmetries coincide only on shell, that is, when the 
equations of motion hold, i.e. in the semiclassical sector of the path integral (critical points of the action). 
But that is hardly surprising. The quantum effects, that is, the fluctuations and higher correlations will 
receive corrections coming from the measure factor and one completely misses them if one postulates the 
naive covariant measure [dq] exp(iS). Notice also that the path integral remembers the gauge fixing condi- 
tion explicitly through the dependence of the wave functions, as well as the exponential of the current, on 
Q A only, rather than all of q a , which are adapted to G. 

4 Operator Constraint Quantisation Path Integral 

As already mentioned, in the presence of second class constraints, operator constraint quantisation is in 
general impossible if one does not pass to local Darboux coordinates with respect to the Dirac bracket 
because otherwise one does not find representations of the canonical commutation relations. Hence we 
assume that we have passed from the global conjugate pairs {q a ,p a ) to local conjugate pairs (z%) = 
(x"",^), (<^,7r^), (Q a ,Pa) in terms of which the second and first class constraints respectively can be 
reformulated as 5s = O z E = and F^ = O F'^ = ir^ + h^((f>; Q , P) = 0. The F'^ are Abelian 
{F'Fl} = and first class {F'[a,zj:} = while {y a ,x a '} = 5%.' and thus the Dirac bracket on functions of 
(j), ir, Q, P reduces to the Poisson bracket. 

On the assumption that the constraints F' can be quantised without anomalies on the kinematical 
Hilbert space Wkin = L,2(dQd§), that is, [F^Fl] = we define a rigging map heuristically as (one has to 
be careful with domains and ranges and 77 should be defined as anti - linear map, however we do not need 
to enter the discussion of these niceties here, see |17| for further information) 

V ■ n kin -> H phys ; V -> foWK&Q) = J WR2ir)} [e^n^Q) (4.1) 

In the case at hand we can easily bring (|4.ip into a form from which it is obvious that it solves the constraints 
F'fj, = 0. First of all we notice that (14. ip can be formally written as 

»7W=IPW ( 4 - 2 ) 

where the order of the 5 distributions is irrelevant due to the Abelianess of the constraints. This is not the case 
for the Fp, which is why the heuristic projector defined in [39] does not solve the constraints. This is already 
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a hint that (|4.ip indeed solves the constraints. To actually prove it we notice that e ll3 ' J ' 7Tf 'ip((p) = ip(<fi — (3) 
and so (we suppress the Q argument in what follows) 



fo(VOK0) = J W V((3) V# - P) = J [d/ty(27r)] V(</> - P) ^(P) (4.3) 

where 

V((3) = e^l^+^W] e" i/3Al7r " (4.4) 

Using 

e^*" ti^e-^* = h' v {<t> - (3) (4.5) 

we derive 

Denoting (3{t) := Pi + t(/3 2 - Pi) it follows 

dt.i 

dt 



V(P 2 ) - V{Pi) = I dh 4-V(P(ti)) = i I dh V(P(h)) $»(0)h'M - Pih)) (4.7) 



where we noticed that $(t) = $(0) = P2 — Pi =const. Equation (|4.7p can be iterated into a Dyson series: 
We need 

PtAh) = Pi + t 2 {P(h) - Pi) = Pi + tMp2 - Pi) = Pihh) (4.8) 



so that 



V(p(h)) - V(P!) = f dt 2 ^-v(p tl (t 2 )) = i f dt 2 V(p(ht 2 )) hp^tiM - p(ht 2 )) 

= i [ 1 dt 2 V(p(t 2 )) P»{Q)h'{4> - P(t 2 )) (4.9) 







Accordingly we obtain for any p 2 , P 



V{px)~ l V(P 2 ) = T, exp(i f dt [p 2 - p x f h'{(j> -fa- t(p 2 - Pi))) (4.10) 

Jo 

where the path ordering symbol T\ orders the earliest time to the left. Choosing /3 X = <f>, P 2 = 4> ~ P we find 

V(<j> -P)= V{<t>) Tt expH I 1 ^ P» h'^tp)) =: V(cp)U(P) (4.11) 

Jo 

where, using again (|4.10p with 2 = (p,Px = and noticing from the definition (14, 4p that V(0) = 1 



V(</>) = T t exp(i f dt p h! - t))) = T r exp(i f dt <p k ' (<frt)) (4.12) 
Jo ' Jo 



where we have performed the change of variables 1 1— ► 1 — t which switches 7} to T r which orders the earliest 
time to the right. For later use we notice the identity 

U(P) = V(P)- 1 = V(P)i (4.13) 

which establishes unitarity of V(P) (as an operator on L 2 (dQ)) and can easist be demonstrated by writing 
V(P) in the form 

V(P) = lim e^WVfl e w^W^p) e £/*%(£/J) ( 4 14) 

N—>oo 
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The point of these manipulations is that we can now write 

i v m(^Q) = v^)wm(Q), W(mQ)= I w/q^vw^^q) (4.15) 



Obviously f]'(ip) no longer depends on (f> so that the rigging map essentially produces functions whose <j) 
dependence is restricted to be of the form V(4>)^>(Q) for suitable ^ G L2(dQ). In order to show that such 
functions really solve = it is very crucial that [F^F^] = 0, otherwise this doe not hold. Essentially, 
the proof boils down to showing (we again suppress the Q dependence) 

[A-- t h>U)}V{4>) = V (4.16) 



This almost looks like a parallel transport equation on (j) space with respect to a one form — i/i^(</>) with 
values in a Lie algebra of (anti - self adjoint) operators on L<2{dQ) defined on a common dense and invariant 
domain and V(4>) looks like its holonomy along the path 1 1— ► t(f>. The difference with the parallel equation is 
of course that the latter is an ODE while (|4.16p is a system of PDE's so that the issue of integrability arises 
and so the following theorem is not trivial (notice tha we do not need to assume d^h' v ] = 0). Its validity 
rests on the fact that 

[F;, F' v ] = -{2d [fM (-ih' u] + [(-*/£), (-<)]} = (4.17) 
i.e. that the curvature of the connection —ih' vanishes. 
Theorem 4.1. 

Equation H^Td ) holds pointwise in cf> space on a dense set of analytic vector£^\ for the operator h'{(j)) = 
</>%(</>). 

Proof. Let V ((f>) := 1 for N G N and for N > 

N /"I I-tn-l 

V N (4>) = 1 + V i n / dh (TKiiht) - / dt n ^ h! Vn {t n 4>) (4.18) 

n=l J JO 

Clearly limjy^oo Vn(4>) = V(4>) converges on analytic vectors for h'{4>) '■= ^h'^cfi). We define for N > the 
remainder 

/"l ftn-1 ftN-1 

/ dh^ti^trf) ... j 

n=l 



Rn{4>) = -i N ' 1 Yl I d *i ^Kxih^) .» / dt n t n ^{[h'^, ^]}M) - J dt N <jr« ti VN {t N cf>) 

(4.19) 

and prove by induction for N > that 

d„V N ((l>) = itipW Vn-M + Rn{4>) (4-20) 

By definition, on analytic vectors of h'(<j>) the norm of the remainder converges (pointwise in <fi space) to 
zero (it is of order 1/[(N — 1)!]) so that once that (|4.20p is established, the proof is complete. 
For N = 1 we obtain 

a M vi($ = i I dtih'^tft + tcpidXKty)] 

Jo 

dt [j f [t h^(t ( p)]+2t^(d [ X ] )(m 

itinW- f 1 dtitp {[h'^h'Mtm 

Jo 

itiJflVofr) + (4.21) 



3 A vector ip is called analytic for an operator A if ||j4. n ^)|| < oo for all n and X^ n =o ^" ll J ^™ 7 /'ll/( n ') < 00 f° r some t > 0. 
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where in the third step we used 

[i^,ifl = & 2id [ll h' u] + W,,K]=0 (4.22) 
Assuming that (|4.20j) holds up to N we compute 

d„(V N+1 - V N ) = i N+l d^ / dh <Th'(ti<l>) ... / dt N+1 <j> v ^h'(t N+1 <i>) 



= i N+1 



Jo 

iV+1 .1 /.*« 



dtx^h'^itxcj)) ... dtnih'^tn^+tn^id^Jitntf))]... dt N+1 0"*+* h' (tjv+K 

n=1 Jo Jo Jo 



JV+1 „1 



^ +1 E /„ dtl <t> VX KM<t>) - /" 1 dtn[^(*n^M)) + 2*n^"(W n ])(W)] 



(i 



n=l " u 
iV+1 „i 



•JV- 



n 



x/ dt n+1 ^h[ ' (t n+1 cf>) ... / rftjv+i 4> VN+1 K ,(^+10) 
Jo JO 

-i N J2 dh^K^) ... dtnt^HK^Kj)^)]... dt N+1 ^+>K N+i (t N+ ^ 

n=l ^° ^° "' 

Eat +1 + ^(^[M^) - Viv-ifa)] 

„= 2 "'o Jo 



dt n+1 ct> u ^h' At n+1 <i>) ... / diiv+i 4> VN+1 h' Attf+u 



II 



/ dh^h'^tKf)) ... ^ 2 dt n -x4f"^h' Vn 1 {tn-l4>) r 1 dt n t n h'^t n <t>)ct> V -h' Vn {t n 4>) X 

„ =1 Jo Jo JO 

x/ dt n+2 4> u ^h'At n+2 4>)... dt N+1 ^ N+1 h'At N+1 <l>) 







Jo 

i?Ar +1 + ^(0)^(0) - v N -A<t>)\ 

N pi rt n -i rtN-i 

/ dh^h'^txcj)) ... dtn0""^M)*n^M)) - / dt N ^ h' VN {t N <t>) 

n=1 Jo Jo Jo 



+* JV - 1 V / dh^K^) ... / 

n=1 Jo Jo 



dt n t n h^{t n (f>) <j> Vn ti Vn {t n 4>) ... / dt N <p» N ti VN {t 



= RN+x+ih'^VN^-VN-A^-RN (4.23) 

In the fourth step we have separated two contributions and the second is easily recognized as the definition 
of -Rjv+i. The non trivial step was the fifth one where we performed an integration by parts in the first 
contribution which produces two sums. We have set to = 1 i n the first sum and in the second in the last 
term the integral over etc. is just unity. In the sixth step we have relabelled in the first sum in the 

n-th term t n+ i — > t n , ..,£tv+i — ► tjy and then reset the summation range to n = 1,..,N. In the second sum 
in the n-th term we have relabelled t n +2 — ► tn+i> --An+i ijv which combines the two sums to —Rn- 
Thus, by assumption (|4.20p 

d„V N+1 = [d^V N - itipVN-i ~ Rn] + ih'^V N + R N+1 = ih'^V N + R N+l (4.24) 

□ 
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Having shown that the rigging map is well defined and produces solutions to the constraints F'^ we 
can compute the physical inner product between states r](ip) defined by (we drop the factors l/(2ir) as the 
physical inner product is defined only up to a scale) 

< ^WO^WO >ph ys := < ^^(ip') > kin 

= j [d<f>] j [dQ]W^Q)V(<f>)W^')](Q) 
= J [d<t>] <^(<t>,.),V(<t>) V '(TP')> L2idQ) 

= < J [^^^-WM^w) 

= < r/(V>W(V0 >L 2 (dQ) (4.25) 

where Vj'{ip) was defined in (|4.15p . This calculation demonstrates that the physical Hilbert space can be 
identified with the Hilbert space 7i Te <± '■= L2{dQ) which we obtained also in the reduced phase space approach. 
The identification is established by 

W : H rcd -> H phys W(Q) i * V(<f>)V(Q) (4.26) 

Wphys can also be recognised as the (closure of the) set of equivalence classes of vectors in 7ikin where ifi ~ ifi' 
iff r]'(ip) = rj'i^') are the same L>2{dQ) functions. Notice that rf is not a projector, [?/] 2 is ill defined. 

It is worthy pointing out the importance of the knowledge of the map (|4.26j) : Often one only knows a 
path integral expression for < r]{ip),r]{ip') > p hys m terms of the boundary states ip,ip' which, however, lack 
any physical interpretation, they are not gauge invariant. The vectors r]{ip) are gauge invariant, however, 
the path integral expression which we will also derive below is not in terms of r](^) but in terms of ip, rf(tp) 
is often not not known explicitly. In the case considered here, r](ip) is known explicitly: Neglecting about 
the details of the domains of the maps we have ^("Hkin) = W(TC re(1 ) and since W just operates by a unitary 
operator with a specific <p dependence, all the non trivial physical information is contained in TC re d- 

To make the link with the path integral formulation now does not require much further work. For any 
G TC re d pick ip,ip' G Wkin with r)'(ip) = ^, rj'ty') = Vl/'. Any such ^ is generated from the cyclic 
vacuum vector (a ground state vector under the time evolution, i.e. a stationary vector under H Te( ±(t) for 
some fixed value of t; in the case of a conservative system, the choice of this t is irrelevant) by operating 
with (limits of) polynomials / of the operators Q A . On the other hand, from the point of view of Ti p h ys 
the operators Q A are ill defined because they are not gauge invariant, or in other words Q A r]{ip) is not 
annihilated by the F'. The following operators, however, preserve Ti. p h ys 

Q A = [exp(i/3^F;)Q A exp(-i^FV)]/3=0 (4.27) 

which is the quantisation of the corresponding classical formula (|2.23p upon replacing {i^,Q" 4 }( n )) by 
[F',Q A ]( n -)/i n . To show that [F^Q^ = we notice that since f^Q" 4 ] = we have with the definition of 

v(p) 63D 

Q A = [expWFU e ~ i/3 ^ Q A e i ^*» exp(-i/3' i F / / u)]/3=0 

= Vi^V^y 1 (4.28) 

Notice that Q A is self-adjoint on H p h ys if Q A is on H re d- Since any physical state is of the form V(4>)^(Q) 
it is obvious that Q A preserves W p h ys by theorem 14.11 We conclude 

< * , *' > red =< n, /(Q)t f(Q)n > red =< wv, w& > phys =< wn, f(Q?f(Q)wn > phys (4.29) 
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We see that the physical scalar product can be directly related to the reduced Hilbert space inner product. 
Now we just need to relate the latter to the n-point functions already derived in the previous section. But 
this is easy: Evidently (|4,29p is a finite linear combination of monomials of the form 

< fl,Q Al ... Q An n > rcd (4.30) 

which is the coincidence limit of an n-point function 

lim < ft,Q Al (*i) ... Q An (t n )n > red (4.31) 
ti,..,t n —*t 

for arbitrary t. In interacting Wightman QFT's it is expected that such equal time correlators are too 
singular [IT]. On the other hand, if the theory can be canonically quantised at all then such limits must 
exist as otherwise the notion of equal time commutation relations is meaningless and therefore presumably 
violates at least one of the Wightman axioms, e.g. the uniqueness of the vacuum. In any case, we derived 
a path integral formula for the right hand side of (|4.31|) in terms of a path integral for the generating 
functional. 

There is also a more direct derivation for a path integral formula for < ry (-(/>), 77 (■;//) > p hys f° r which, 
however, the relation to the reduced phase space path integral is less clear. On the other hand that 
alternative derivation makes the connection to the Master constraint path integral clearer. We will thus 
display it here for completeness. We start from the definition of the rigging map (|4.ip . choose some arbitrary 
but fixed reference vector f^o and normalise the physical inner product by asking that the norm of rj(ip) be 
unity. Thus we have to divide (|4.ip by a constant up to which the inner product is anyway undetermined 
and obtain 

f[d(3] <ii,e^ F ^'> kh 
f[d0] <n ,e^ F ^Q > ki _ 

Notice that (|4.32p is not a path integral over 0, it is just an integral at fixed time of the Lagrange multipliers 
(3^ . In order to introduce a path integral of Lagrange multipliers we introduce an arbitrary time parameter 
T which we will eventually send to oo and multiply both numerator and denominator of (|4.32p by the infinite 
constant 

C = I [DX] Y[6(f dtX»(t)) (4.33) 

which is a path integral over paths t t— > X(t), t E [—T,T]. By shifting the integration variable X(t) = 
X'(t) — 7^f/3 for any constant path /3/(2T) we see that C can also be written 

C = ! [DX] Yl S( f dtX^{t) - (3") (4.34) 

where (3 is arbitrary. Inserting this into (|4.32|) and interchanging the [DX], [d/3] integrals we obtain 

/ [d(3] < ^e^V >kin [/ [DX] EL SU-t dtX»{t)-{3»)] 



. (is J i u h> ^ W,^ >-W ^kin (aqo\ 
< VW,VW ) > P hys= rrin , — iMF'n ( 4 ' 32 ) 



< r)(ip),r)(ip') >p h: 



/ [d/3] < n ,e^J2o > kin [J [DX] n M 5(f T T dt\»(t) - ^)] 

j [dx] < ^^-t > kiD [j m s(j T T dtm^-jn 

J [DX] < n ,e^-T <^«fi > kin [J [d(3] n„ 5(f F T dtX»(t)-(3»)\ 
J[DX] <^,e^-T dt ^ F '^i >]dn 

j [dx] <n ,jV-T^(mn 0>kiI1 



(4.35) 



By writing 

[f dtX»(t)]F'^ = ]imN^oo± jr X»{nT/N)Fl (4.36) 

~~ ' n=—N 



rT i AT-1 
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we finally obtain using the usual skeletonisation techniques 



, „ , „, / [DQ DP D(f> Dtt DX] V(Qt, <h) *P'(Q-T, 4>-t) e< dt[PA ^ +7r " W^*)] 
< rf(ip),rf(ip ) >= -y — : 

/ [DQ DP D4> Dtt DX] tp(Qr, <h) ViQ-T, 4-t) e* Lt dt[PA Q +7r " P-W^QM,*)] 

(4.37) 

Notice that all canonical coordinates and Lagrange multipliers are integrated over paths in time in the 
interval [— T, T] and that the operator F' has been replaced by the classical function in (|4.36p . In this 
expression the parameter T is arbitrary and we can take T — > oo. 

In order to invoke the gauge fixing conditions = T**(i) — <^ we will make use of the Fadeev - Popov 
procedure: Let cty(t) = exp(7 M (t){F^(t), .}) where F'(t) = F'(u(t)) is the constraint on the copy of the 
phase space at time t and u(t) = (Q(t), P(t), cj>(t), ir(t)). Then a l{t) {G^(t)) = G»(t) +^{t) so that in this 
case trivially 

"[Dj] J] 5(a l{t) (G»(t)) = l (4.38) 

We multiply both numerator and denominator of (|4.37p by this unity. Assuming limy^oo j(±T) = the 
kinetic term in the exponential of (|4.37p is invariant (being a symplectic potential), F'^ is invariant due to 
the Abelianess and the Liouville measure at time t is invariant under the canonical transformations a 7 (i). 
Thus after a change of variables from u — > a 7 (u) and since [a 7 (G)](it) = G{a 1 {u)) nothing depends on 7 
anymore and the integral over D'j can be dropped. We obtain 



/ 



/5 T dt[p A 



< ,«■),,»<) >= - lDQ D P D * PT| m**!* S ™ . ■ (4.39) 

S! «2t, <h) "oW-I-, 0-t) i[Gl e''-' "' [ ^ °' 4+ *» *"] 

where we have also integrated over A. 

Finally, in order to invoke the second class constraints in the form z-^ = (x a ,y a ) = we simply insert a 
5 distribution 5[z] and integrate over z. This yields 

< M) M') >= 1 [DQ DP D<P — Dy] 0T) <A ~ T) S[F ' ] 6[G] 6[Z] dt[PA QA+ ^ ^ +V ^ ] 

/l " /l ' / [DQ DP Dcf> Dtt Dx Dy] O O (Q T ,0 T ) fi (Q-r, 4-t) S[F'] S[G] 5[z] e^-r dt ^ A 

(4.40) 

Next we observe that det({F', G}), det({z, z}) are constant in the system of coordinates chosen so we can 
multiply numerator and denominator of (14.401) by these constants. As established in section O the expression 

5[F'] 5[G] 5[z] I det({F', G})\ ^det({z,z} (4.41) 

is invariant under any mapping (F',G,z) 1— > (F,G',S) as long as both triples reduce to the same gauge 
cut of the same constraint surface. We may therefore restore the original first and second class constraints 
F, S while keeping G = G' provided we keep the determinant factors in (|4.41jl . Finally we can restore the 
original system of coordinates q a ,p a which arise from (q' a ,p' a ) '■= (Q a ,Pa), {<P^iT^u)i (x a ,ya) by a canonical 
transformation a because the symplectic potential in the exponent of (|4.40p as well as the Liouville measure 
remain invariant and the the Poisson brackets are simply expressed in the new coordinates, e.g. 

{S,S}(q',p') = {S,S}(a(q,p))={Soa,Soa}(q,p) (4.42) 

(by S we denote the original S expressed in whatever canonical coordinates). Accordingly 

/ [Dq Dp] 4>{Qt, 4>t) WQ-tA-t) S[F] 5[G] S[S] |det({F,G})| y/det({S,S}) 



[ [Dq Dp] n (Q T ,cl>T) fi o (Q-T,0-T) S[F] S[G] 5[S] \det({F,G})\ y/dek({S,S}) e^-* dtpa ' a 

(4.43) 
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Notice that due to the gauge fixing condition G{t) = r(i) — <f>{t) the integral over <p is anyway concentrated 
at the fixed path r(t) so that it is allowed to assume that ip , t/j' , f2 are actually independent of <j). In this 
sense the final result (|4.43j) precisely agrees with (j3.70p with the understanding that ifj, ip' in (|4.43|) can be 
generated from the generating functional (|3.70p by suitable functional differentiation with respect to the 
current j at j = at coincident points of time ±T in the limit T — > oo. 

5 Master Constraint Path Integral 

The Master Constraint Programme (MCP) was originally designed precisely in order to be able to cope with 
gauge systems whose classical first class constraint algebra involves structure functions [16] and for which 
therefore group averaging techniques do not work. It is true that locally the first class constraints F can 
be replaced by equivalent ones whose algebra is Abelian and we have made heavy use of that fact in the 
two previous sections. However, for the case of interest, namely General Relativity, in vacuum the Abelian 
constraints are rather non local on the spatial manifold, algebraically difficult to deal with and not explicitly 
known even classically [40]. Even with standard matter this is true. It is for this reason that in [10^ [33] 
non standard matter (Brown - Kuchaf Dust [U]) was used in order to achieve the Abelianisation in a local 
form and such that the resulting expressions remain practically managable. The MCP does not rely on 
Abelianisation and thus is both more global (on phase space) in character and does not require any special 
type of matter. In principle it does not even require that the constraints are quantised without anomalies 
and even second class constraints can be treated by the MCP [16] . Since the Master Constraint is a weighted 
sum of squares of the first class constraints, we expect that its kernel is empty when the constraints are not 
quantised without anomalies. In that case one could consider the Hilbert "subspace" corresponding to the 
lowest "eigenvalue" as the suitable substitute for the anomaly free situation. See [16] for further discussion. 
In that sense the MCP may be considered as a much more flexible approach to constrained systems with 
structure functions. 

While for a wide range of models the MCP has been tested versus the more traditional operator constraint 
method [16], its equivalence with the latter is so far lacking. On the one hand, the equivalence seems to 
be obvious since both the Master constraint and the individual constraints are supposed to define the same 
(common) kernel. On the other hand, the equivalence is rather not obvious because the formulae for defining 
the physical inner product or equivalently the rigging map are totally different. For the individual constraints 
in Abelianised form the rigging map is defined in (|4.ip while for the MCP it is heuristically defined bJ^l 



The symmetric (possibly operator valued) matrix K should be so chosen such that M is positive and such 
that it arises from a classically positive definite matrix valued function on phase space. There is great 
flexibility in the choice of K and while all (sufficiently differentiable) positive definite classical matrices are 
equivalent, in quantum theory this flexibility must be exploited in order to arrive at well defined master 
Constraint Operators [16] . Normally we require that F^ is quantised as a self adjoint operator but in the 
case of structure functions this must be relaxed [16] which is why we included the adjoint in (|5.2p . 

The task of the present section is to connect with the results of the previous two sections. Those 
sections made use of the Abelian constraints F' and we will therefore use those in order to build our Master 
Constraint. We assume as in sections [3] and that F'^ is self - adjoint since the F' are supposed to be 

19 Again there are subtle domain issues which we neglect here and moreover one should switch to to a direct integral repre- 
sentation of 7ikin subordinate to M; see |16| for details. 




(5.1) 



where the Master constraint is defined by 




(5.2) 



36 



quantised without anomalies. As in the previous section we choose a reference vector and define the 
Master Constraint physical inner product by 



M J, dt <V,e«"V'> kl , 



To see that (14. ip and (15. 3p formally coincide, recall [1] that for any self - adoint operator A on a (separablJ^l) 
Hilbert space there exists a unitary transformation (generalised Fourier transform) 



U : u -» W® := / d)«(A) W®; V ^ ty>(A)) Aespoc(A) (5.4) 

«/ spec(/l) 

from to a direct integral of Hilbert spaces W® (possibly with different dimensions for each A but in a 
maeasurable way, hence more general than a Hilbert bundle) with respect to a probability measure /i on the 
spectrum spec(^4) of A. Here V>(A) € W®. The correspondence between the inner products is 

< W >h=< 4>,fi > H ®--= / < ^(A),^'(A) > Hffi (5.5) 



spec(yt) 



The point of this spectral decomposition is that [£/.A?7 _1 ^;](A) = A^(A), i.e. A acts by multipliction by A 
on 7if If (the spectral projections of) two self - adjoint operators A, B commute then UBU^ 1 preserves 
H'f and we may apply the just quoted theorem which then tells us that there exists a joint probability 
measure cI(i(\a, Xb) on the joint spectum spec({A, B}) = spec(vl) x spec(-B) of A, B and a representation of 
H as a direct integral of Hilbert spaces 7i®^ Afl on which A, B respectively act by multiplication by Xa, Xb 
respectively. 

Iterating like that we obtain the statement that for a (countable) family of mutually commuting self - 
adjoint operators F' there exists a unitary operator U from Wkin to 7i® which is the direct integral with 
respect to a measure [i on the joint spectrum of the of Hilbert spaces W® Am}m on which UF^U^ 1 acts by 
multiplication by A^. This is the key to link (|5,3p and (|4.ip . Namely we formally obtain for (|4.ip 



phys 



f[d/3] <fi ,e i/3M ^ Q 0>kiD 

L P ec( { F> }) MW) <^({A»,^({A})> M e } [J 

Jspcc( { F'}) MW) <n"o({A}),n ({A})> w e } [fM^i 

/s P ec({F'}) ^M({A}) <^({A}),^({A})> w e } 5({A}) 
/spcc( { F'}) ^({A» <fio({A}),n ({A})> H e } «5({A}) 
p({0}) <^({0}),#({0})> < 

{0} 

p({0}) <fi ({0}),Oo({0})> w e 

{0} 

<^({0}),^({0}) > w e } 



<O ({0}),Q ({0}) > Hf 

{0} 



(5.6) 



where formally dfx({X}) =: p({A})[<iA]. Notice that p({A}) can have distributional contributions if the 
spectrum has a pure point part, see [T6J ES]- Of course there are measure theoretic issues such as: if lies 
in the continuous spectrum of some F^ then {0} has fi measure zero and is not well defined without 



20 In LQG the Hilbert space is not separable but the operator M preserves the separable subspaces into which the Hilbert 
space decomposes 
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further assumptions spelled out in [16] . For the purposes of this paper we take a formal attitude and simply 
let the formal cancellation of the p({0}) in numerator and denominator of (|5.6p take place as indicated. For 
a more careful definition see |25| . 



On the other hand we have 



J dt <V>,e* tA V> kin 

/ dt < n ,e itM n > kin 

JspecUF'}) MW) <^({A}),^({A})> H e } [/ dte U ^ K ^ X ^] 
L P ec( { F' } ) MW) <no({A}),fio({A})> w e } [J dte u ^ K ^ X »] 

Jspec( { F' } ) MW) <^o({A}),^ ({A})> w e } SCE^K^XM 
p({0}) J({0}) Vol(5) < ^({0}),^'({0}) > H? 

{0} 

P(W) J({0}) Vol(S) < fi o ({0}),fi ({0}) > H e o} 
<^({0}),^'({0}) > w ® } 



<fi ({0}),fio({0}) > w e 

{0} 



(5.7) 



where J({X}) is the Jacobian that arises by switching from {A} to polar coordinates adapted to the radius 
squared r 2 := E^ v K^X^Xu. Of course we have assumed that is just a complex valued positive definite 
matrix. Vol(S') is the volume of the corresponding sphere. Fo countably many F'^ the volume of the infinite 
dimensional sphere vanishes as well as the Jacobian at zero. To justify (15. 7|) less formally one has to take a 
limit as the number N of F^ approaches infinity so that Vol(S' Ar_1 ) is finite and one also has to regularise 
5{M) by 8{M — e 2 ) and take e — > as to make J(e) finite. See [25] for the details and also (15.10P below for 
a sketch. 

Hence (j5.6[) and (|5.7p agree with each other modulo formal manipulations and thus give rise to the same 
path integral formulation. Our method of "proof" above used spectral theory. We will now provide a more 
direct (but also formal) "proof" using only path integral techniques. The idea is the same as at the end of 
section [5] and was already sketched in [16] . First of all we use the same technique as used between (I4.32p 
and (I4.37P in order to write (15. 3D as 



<wW-1mW > M 



phys" 



J[dx] < ^a {Lt dt mw^ >k . 



J[DQ DP D<j> Dtt DX] ^(QtAt) ^'{Q-tA-t) dtKP A Q^Wt)-\(t)Mm)W)Mt)rtt)] 
J[DQ DP Dcf> Dtt DX] ^{QtAt) fi (Q-T, 4>-t) e if ~ T *K^O A -H»r»(*)-A(*)M«(*),P(t) I *(t), 1 r(t)] 

\T] 



f[DQ DP D<p Dtt] ^(Q T ,<h) ^'(Q-tA-t) [Ute[-T,T] <*((Af(i))] e ^-r dt[{P A Q A +^){t)] 

if_ T dt[(P A QA +7T ^)(t)] 



f[DQ DP D(p Dtt] O o (Qt,0t) ^{Q-tA-t) [lL G [-T,T] S((M(t))} e 



(5.8) 



where T is again an arbitrary parameter which we take to oo eventually. If in (|4.37p we perform the 
integral over A then the only difference between (|4.37|) and (|5.8|) is that instead of 5[F'] the distribution 
S[M] appears in both numerator and denominator. But clearly the two distributions have the same support 
7r = —h'(4>,Q,P). Let us therefore explicitly do the integral in both (|4.37p and (|5.8p and compare the 
results. It suffices to do this at fixed t because both 5 distributions factorise over [— T, T]. We consider S(M) 
as the limit iV — ► oo, e — > of 

S N>e (M):=S( J2 K^F^Fl-e 2 ) (5.9) 
H,v<N 
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Let / = /[tt] be any functional of tt^, /i = 1,..,N. The N x N submatrix = K^ u \ /i, v < N is also 
positive definite on the corresponding vector subspace. Hence its square root and inverse is well defined. 
Thus, by shifting the integration variable and switching to radial r and polar coordinates ip respectively we 
obtain with the unit vector x^/r = n^tp) 



d N 7T 5 N>e (M) /(tt) 

f d N x 5(x T K N x - e 2 ) f(-ti + x) 
Jm N 

- 1 f d N x 5{x T x - e 2 ) f(-ti + K h 
^det(K N ) Jrn 



1/2, 



1 



r^" 1 dr 5{r 2 - e 2 ) / dV61(y>) f(-ti + K~ 1/2 rn(ip)) 



y/det(K N ) Jr + Js" 



JV-2 



2y/det(K N ) JsN-i 
In the limit e — * this approaches 



dV61(y>) /(-/t' + K N 1/2 en(ip)) (5.10) 



2y / det(^iv) 
and in that sense we may write 

JV-2 



Af-2 

Vol^- 1 ) /(-^ (5.11) 



6 N , e (M) = r—— VoK^" 1 ) 5n(F'), S N (F') = ]J 5(F;) (5.12) 
V det (-^iv) 

Since i^Ar is a phase space independent constant, when inserting (|5.12p into (|5.8p . the prefactor cancels in 
both numerator and denominator and we arrive at (|4.37p in the limit e — » and N ^ oo. 



6 Conclusions and Outlook 

The three tasks accomplished in the present paper are: 

1. We have demonstrated that within the limits of the formal nature of the manipulations that are usually 
employed when dealing with path integrals, three canonical quantisation methods, namely the reduced 
phase space -, the operator constraint and the Master Constraint quantisation all lead to the same 
path integral formulation for the physical inner product. In order that rigging map techniques can be 
employed to the operator constraint approach, in the case of structure functions one has to pass to 
Abelianised constrants. 

2. The resulting path integral can be written in terms of the classical Lagrangian from which the classical 
theory descends. However, the correct measure to be used is not the naive Lebesgue measure on path 
space, rather this measure must be corrected by factors that depend on the first and second class 
constraints as well as the gauge fixing condition. 

3. The gauge fixing condition is in one to one correspondence with the choice and interpretation of a 
convenient choice of an algebra of physical observables and a physical Hamiltonian. It is possible to 
do without gauge fixing conditions provided one finds alternative methods to construct an algebra 
of Dirac observables. However, the resulting algebra is almost surely algebraically more complicated, 
more difficult to quantise, lacks an a priori physical interpretation and is not equipped with a preferred 
physical time evolution. In particular, if one wants to talk about the scattering matrix between physical 
states, the dependence on the gauge fixing is unavoidable because it determines the physical time 
evolution of the chosen "basis" of gauge invariant operators. 
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As we have already stated in the introduction, certainly not all the results and techniques derived and used 
in the present paper are new, bits and pieces of it are already in the literature. However, we believe we 
have assembeled the material in a new and fruitful way in order to better understand the relations between 
the four quantisation methods discussed in this paper. Also we think that the mathematical and physical 
influence of the gauge fixing condition hasbeen described in this paper from a new angle. 

As we have seen explicitly, both methods of proof in section [5] actually relied on the fact that the matrix 
K is a constant function on phase space. However, this is not the case for the concrete Master Constraint 
for General Relativity studied in [16J. Namely, there one considered an expression of the form 

r C 2 

M = d 3 x —== (6.1) 
J a v / det(g) 

where C is the Hamiltonian constraint and q is the intrinsic three metric of the hypersurface a. The "matrix" 
K(x,y) = 6(x,y)/ y' det(q)(x) is chosen here in order to make (|6.ip invariant under spatial diffeomorphisms 
and is clearly a non - trivial function on phase space. In view of the analysis of the previous section, rather 
than the Hamiltonian constraint in its original form C, in the presence of the dust matter one would choose 
it in the locally equivalent form C'(x) = ir(x) + h'(q(x), P(x)) where q,P are the gravitational degrees 
of freedom and it is one of the dust momenta. Notice that for this type of matter h! does not depend 
on the dust configuration fields (f> and therefore dust deparametrises the system and leads to a conserved 
physical Hamiltonian. However, also C is a scalar density and thus to make the corresponding Master 
constraint spatially diffeomorphism invariant, one would again need a phase space dependent matrix of the 
type considered above. Thus it appears as if the analysis of the present section does not apply to GR. 

However, this is not the case: Namely, the dust offers the possibility to completely abelianise the full 
constraint algebra including spatial diffeomorphisms. Thus in contrast to the usual situation in which the 
spatial diffeomorphisms form a subalgebra of the constraint algebra but not an ideal, it is possible to com- 
pletely solve the spatial diffeomorphism constraint before solving the Hamiltonian constraint. In particular 
it is possible to perform a canonical transformation to coordinates such that C only depends on spatially 
diffeomorphism invariant fields |10j . It is therefore no longer necessary to choose a density weight minus 
one matrix K. We can simply take an orthonormal basis of L2(cr, d 3 x) and consider the F' :=< b^, C >. 
Then one chooses any phase space independent matrix K^ v subject to certain fall off conditions in index 
space (typically K should be trace class |16|). The fact that C has density weight one ensures that C can 
be quantised on the unique [12] LQG Hilbert space [13] as was shown explicitly in [33]. That quantisation, 
however, is most probably too naive in order guarantee anomaly freeness and must be improved. Yet, since 
the anomaly is an h correction to the classical result, the relation between the MCP (which also works in 
the anomalous case) and the path integral formulation derived in the previous section, remains correct in 
the semiclassical limit. An alternative to working with C already reduced with respect to the spatial dif- 
feomorphism constraint is to keep the unreduced C and the unreduced Abelianised spatial diffeomorphism 
constraints Cj |10] . The caveat in LQG to quantising the classical generator of spatial diffeomorphisms 
which arises due to strong discontinuity of the one paranmeter unitary subgroups of spatial diffeomorphisms 
on the LQG Hilbert space is circumvented because Cj is not a density one covector but a density one scalar 
and thus can be quantised on the LQG Hilbert space [33], albeit it is difficult, similar to C, to achieve 
anomaly freeness. 

This paper has been the starting point for further analysis. In [21] we have computed the correct mea- 
sure for the Hoist action and have checked explicitly that it is consistent with the analysis of [20] for the 
Plebanski action. In |25j the relation between the Master Constraint Programme and the operator con- 
straint programme for Abelian and anomaly free constraints and with phase space independent matrix K 
was analysed with higher mathematical precision at the level of the canonical theory and it is shown that 
under certain technical assumptions the two methods lead to the same result, thus partly removing the 
formal character of the analysis of section [5j Finally, in [26] it was formally checked by using available 
semiclassical techniques |44] that the Master Constraint Programme for General Relativity leads also to 
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the expected path integral formula up to a local measure factor when one considers phase space dependent 
matrices K and non - Abelian constraints. However, the results here are less strong (more formal) than in 
the Abelianised case. 

Many further questions arise from the present paper: 
Since the Master Constraint can in principle also accomodate (sums of squares of) second class constraints 
if one subtracts a suitable normal ordering constant [16], one could ask whether the separate treatment of 
first and second class constraints could be unified and if yes how the corresponding path integral would look 
like. Secondly, in applications to path integral formulations of LQG one should really take the unavoidable 
measure factor derived in [21] and following the general theory summarised here seriously and define a 
corresponding spin foam model. Work is in progress in order to achieve that. Next, due to the measure 
factor the theory lacks manifest spacetime diffeomorphism invariance. On the other hand it should be man- 
ifestly invariant under the gauge transformations generated by the first class constraints which in General 
Relativity corresponds to the Bergmann - Komar "group" [36] (more precisely it is the enveloping algebra 
generated by the secondary first class constraints of GR). The two groups are known to coincide when the 
classical equations of motion hold and this is the reason why the Lagrangian and Hamiltonian descriptions 
are equivalent classically. However, off shell there is no particular relation between these two "groups" and 
it is consistent with the classical theory that the spacetime diffeomophism group is not a symmetry of the 
quantum theory. In [22] it is further analysed in which sense the Bergmann - Komar group is a symmetry of 
the Hamiltonian path integral. It seems that the attempt to construct a spacetime covariant path integral 
of GR has no chance to be derived from a canonical platform which is the only systematic starting point 
that we have and it would be interesting to understand better the implications of this conclusion. In some 
sense it is clear that spacetime diffeomorphism invariance is far from sufficient in order to guarantee that 
one has a correct quantisation of a given classical theory. Many Lagrangians are spacetime diffeomorphism 
covariant (e.g. higher derivative theories) but all of them have different Hamiltonian constraints (even dif- 
ferent numbers of degrees of frredom). The effect of this will show, in particular, in the local measure factor 
that we have exhibited. 
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